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SUMMARY  OF  RESEARCH  ACCOMPLISHMENTS 


This  research  has  developed  a  solution  procedure  for  a  practical,  new  combinatorial 
optimization  problem,  namely,  problem  (*)  of  page  4  of  this  technical  report.  The  solution 
procedure  also  represents  a  new  approach  to  combinatorial  problems. 

Problem  (*)  is  a  model  for  certain  surveillance  maintenance  problems,  and  has  not 
previously  appeared  in  the  literature  of  optimization  or  industrial  engineering.  Therefore,  this 
research  is  bringing  this  important  problem  to  the  attention  of  researchers  and  specialists  in  the 
field  of  optimization  and  industrial  engineering. 

Our  procedure  for  solving  problem  (*)  is  a  heuristic  that  produces  good  quality, 
approximate  solutions.  This  heuristic  is  designed  to  utilize  solution  guesses  from  the  user,  and 
incorporate  new  fast  heuristics  for  the  standard  set-covering  problem.  Hence,  this  procedure  has 
the  potential  to  "evolve"  and  "develop",  as  more  computational  experiments  are  performed,  and 
researchers  become  aware  of  problem  (*). 

This  heuristic  introduces  a  new  approach  to  combinatorial  problems.  Combinations  of 
sensor  mechanisms  are  "lined  up"  as  (binary  coded)  integers,  and  then  a  discrete  line  search  (that 
is,  a  line  search  on  integers  only)  is  performed.  Hence,  a  discrete  optimization  problem  is 
transformed  in  a  way  that  allows  basic  procedures  of  nonlinear  programming  to  be  applied.  This 
should  be  a  useful  approach  for  a  host  of  NP-complete  combinatorial  problems  (the  knapsack 
problem,  for  instance). 

Finally,  the  cubical  subdivision  of  the  set  S  described  in  Appendix  1  of  this  report  should 
have  interesting  implications  for  research  on  surveillance  sensors. 
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SECTION  1 


INTRODUCTION 

This  is  the  final  report  for  research  contract  DASG60-90-C-0142  with  the  United  States 
ArtTiy  Strategic  Defense  Command,  Huntsville,  Alabama.  This  is  also  the  third  of  three  technical 
reports  on  the  research  contract.  The  first  and  second  preliminary  technical  reports  are  included 
in  this  final  report  as  Appendix  I  and  Appendix  II  respectively. 

The  subject  of  this  research  contract  is  the  problem  labeled  (*)  below: 

(*)Given:  S,  a  subset  of  the  3-dimensional  space,  R’;  subsets  Sy,  i=l,...,m,  j=l,...,n,  that  cover 
S  in  the  sense  that  S  is  the  union  of  all  the  S^’s;  two  numbers  tjj  and  c^  that  will  be 
interpreted  as  "time"  and  "cost"  respectively.  Required:  To  determine  an  optimal  collection 
of  the  Sjj’s  that  cover  S,  with  no  more  than  one  Sy  for  each  i,  and  optimality  defined  in 
terms  of  minimizing  max{tij  :  Sy  is  in  the  collection)  or  minimizing  the  total  cost  of  the 
Sy’s  in  the  collection. 

This  set-covering  problem  shall  be  referred  to  as  problem  (*).  Problem  (*)  is  a 
representation  for  a  class  of  surveillance  maintenance  problems  including  the  problem  of  restoring 
surveillance  coverage  for  a  system  of  surveillance  sensor  mechanisms  after  a  sensor  mechanism 
(in  the  system)  fails.  The  set  S  represents  the  space  under  surveillance;  each  subset  Sy  represents 
the  portion  of  S  covered  by  the  i-th  surveillance  sensor  mechanism  functioning  in  its  J-th  optional 
deployment/positioning;  each  Cy  represents  the  cost  of  the  j-th  deployment  of  the  i-th  sensor 
mechanism;  tjj  represents  the  time  required  for  the  i-th  sensor  mechanism  to  be  placed  in  its  j-th 
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optional  deployment.  Selecting  a  minimum-cost  set  of  Sy’s  covering  S  corresponds  to  restoring 
surveillance  coverage  at  minimum  cost,  whereas  selecting  a  minimax-time  set  of  Sy’s  covering 
S  corresponds  to  restoring  surveillance  coverage  within  the  shortest  time  possible. 

For  any  procedure  to  solve  problem  (*),  it  is  necessary  to  be  able  to  determine  whether 
a  given  selection  of  the  Sy’s  define  a  cover  for  S.  This  question  is  the  subject  of  our  first 
preliminary  technical  report.  In  addressing  this  question,  we  find  information  about  the  sizes  of 
potential  intruders  (into  S)  to  be  very  helpful.  This  information  allows  us  to  avoid  putting 
stringent  conditions  on  Sy  or  S.  Section  2  of  this  report  gives  a  summary  of  the  preliminary 
technical  report. 

Also  in  section  2  of  this  report,  our  solution  procedure  for  problem  (*)  is  summarized. 
Appendix  II  contains  a  detailed  description  of  this  procedure.  Partly  because  problem  (*)  is  yet 
to  appear  in  established  literature,  our  procedure  for  solving  problem  (*)  is  not  a  variant  of  any 
existing  procedures.  Section  3  of  this  report  evaluates  certain  computational  features  of  this 
procedure. 

Sections  3  and  4  of  this  report  were  not  included  in  the  preliminary  technical  reports. 
Using  various  test  problems,  section  3  contains  information  on  the  quality  of  solutions,  and 
discusses  some  design  choices.  The  computational  complexity  of  the  procedure  is  already 
estimated  in  Appendix  II  as  o(ISIm^log,  32"'). 

Section  4  describes  three  extensions  of  problem  (*)  as  possibilities  for  further  work. 
These  are  resource-constrained  extensions,  commitment-constrained  extensions,  and  multiple-layer 
coverage  extensions.  Such  extensions  incorporate  a  more  realistic  setting  for  problem  (*). 

To  close  this  introduction,  a  clarification  of  one  terminology  that  will  be  encountered 
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frequently  in  this  report  is  in  order.  Since  it  is  a  heuristic,  our  solution  procedure  is  application- 
oriented,  and  so  is  the  description  of  it.  Hence,  sometimes  problem  (*)  is  referred  to  as  the 
surveillance  maintenance  problem,  even  though  problem  (*)  is  actually  an  optimization  model 
that  has  the  surveillance  maintenance  problem  as  its  motivation,  and  one  of  several  possible 
applications. 
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SECTION  2 


DESCRIPTION  OF  SOLUTION  PROCEDURE 

Our  solution  procedure  has  been  formally  described  in  our  first  and  second  preliminary 
technical  reports,  DASG60/90/1  and  DASG60/90/2,  which  are  included  here  as  Appendix  I  and 
Appendix  II  respectively.  In  this  section,  we  describe  our  motivation  for  the  results  presented 
in  those  preliminary  technical  reports.  Informal  summaries  of  the  results  are  also  given. 

2.1  Determining  Whether  S  is  Covered 

Any  method  for  solving  problem  (*)  must  include  a  procedure  for  checking  whether  the 
set.  S  is  equal  to  the  union  of  any  selection  of  the  Sy’s,  that  is,  whether  S  is  covered  by  a  given 
selection  of  the  Sy’s.  If  the  set  S  were  a  finite  set,  then  it  would  be  a  straightforward  matter  to 
do  this  checking  on  a  digital  computer.  But,  in  problem  (*),  the  set  S  in  not  finite;  it  is  a  3- 
manifold  in  (that  is,  a  full-dimensional  volume  in  R^).  This  task  of  checking  whether  S  is 
covered  is  the  subject  of  our  first  preliminary  technical  report. 

In  the  case  that  S  is  not  a  finite  set,  it  may  be  very  difficult  or  even  impossible  to  check 
whether  S  is  covered.  To  illustrate  this  point,  let  us  consider  two  special  examples.  First, 
suppose  S  and  the  Sy’s  are  rectangles  in  R^.  Figure  2.1  below  portrays  a  possible  situation  that 
may  develop.  This  diagram  suggests  that  the  task  (checking  whether  S  is  covered)  is  not 
impossible  in  this  case,  but  could  involve  a  great  deal  of  record  keeping. 

Now  suppose  S  is  a  rectangle  in  R^,  and  some  of  the  Sy’s  are  intersections  of  circles  and 
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S.  Figure  2.2  below  portrays  a  possible  picture  in  this  case.  This  diagram  suggests  that  it  may 
be  impossible  to  check  whether  S  is  covered,  unless  further  information  is  brought  to  bear  on  the 
case. 


Figure  2.2 


In  problem  (*),  S  corresponds  to  the  space  under  surveillance,  and  each  Sjj  corresponds 
to  the  portion  of  S  covered  by  the  span  of  some  sensor  mechanism  (in  a  specified  deployment). 
Sjj  depends  on  sensor  technology,  whereas  S  is  a  fixed  volume.  Since  sensors  have  various  non¬ 
surveillance  uses  in  the  real  world,  it  may  not  be  realistic  to  impose  restrictions  on  the  Sjj’s  just 
to  make  it  convenient  to  check  whether  S  is  covered. 
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Fortunately,  the  nature  of  surveillance  maintenance  itself  provides  a  useful  information, 
namely  ,  the  size  of  a  potential  intruder.  Suppose  x  is  the  size  of  the  smallest  potential  intruder 
into  S.  S  may  be  "subdivided"  into  a  finite  number  of  pieces  that  fit  together  suitably,  each  piece 
of  size  at  most  x.  Then  the  surveillance  sensor  mechanisms  only  have  to  cover  (all  of)  the 
vertices  (comers)  of  those  pieces.  Since  there  is  a  finite  number  of  such  pieces,  this  essentially 
reduces  S  to  a  finite  set  for  the  purpose  of  comparing  S  to  the  union  any  given  selection  of  the 
Si/s. 

In  the  first  preliminary  technical  report,  alternative  ways  of  subdividing  S  for  this  purpose 
were  considered,  and  a  particular  kind  of  subdivision,  the  cubical  subdivision,  was  found  to  be 
the  most  suitable.  The  cubical  subdivision  allows  an  easy  denumeration  (line-up)  of  the  pieces, 
and  does  not  involve  too  many  vertices.  Other  desirable  properties  of  the  cubical  subdivision 
were  also  described  in  the  technical  report. 

The  literature  of  piecewise  linear  topology  contains  various  types  of  such  a  subdivision 
of  manifolds.  Subdivisions  of  manifolds  are  also  used  extensively  for  homotopy  methods  in 
mathematical  programming  [1][19]  and  finite  element  methods  [15] [18]. 

2.2  Our  Procedure  For  Problem  (♦) 

It  has  been  shown  that,  for  the  purpose  of  solving  problem  (*),  S  may  be  regarded  as  a 
finite  set  (by  virtue  of  results  given  in  Appendix  I).  But,  even  with  S  regarded  as  a  finite  set, 
problem  (*)  remains  ordinarily  harder  than  the  standard  set-covering  problem,  since  problem  (*) 
reduces  to  the  standard  set-covering  only  when  n=l  and  cost  is  the  objective  function.  Therefore, 
one  may  not  rely  on  standard  set-covering  procedures  to  solve  problem  (*).  Our  procedure  for 
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solving  problem  (*)  is  the  subject  of  the  second  preliminary  technical  report  (Appendix  11). 


Problem  (*)  is  an  NP-complete  combinatorial  problem,  as  the  standard  set-covering  is  NP- 
complete.  By  the  theory  of  algorithms  [7]  [14],  this  means  that  it  is  highly  unlikely  that  anyone 
can  ever  find  an  "efficient"  procedure  that  computes  an  exact  solution  for  problem  (*).  As  an 
illustration  of  the  meaning  of  this  NP-completeness  statement,  solving  an  instance  of  problem  (*), 
with  10  sensor  mechanisms  and  5  deployment  options  for  each  sensor  mechanism,  may  take  years 
on  a  fast  digital  computer  doing  one  operation  per  micro- second,  if  an  exact  solution  is  required. 
Hence,  a  realistic  approach  for  real  world  problems  is  to  use  procedures  that  produce  good 
approximate  solutions.  Such  approximate  procedures  are  known  as  heuristic  procedures  or, 
simply,  heuristics. 

Since  they  are  designed  to  produce  approximate  solutions,  rather  than  exact  solutions 
(even  though  they  may  generate  exact  solutions  quite  often),  heuristics  are  usually  designed  to 
meet  quality  specifications  derived  from  particular  applications.  Indeed,  it  is  generally  believed 
that  any  heuristic  must  be  motivated  from  the  needs  of  some  specific  applications. 

Accordingly,  our  procedure  for  solving  problem  (*)  is  a  heuristic  whose  features  are 
motivated  by  the  needs  of  surveillance  maintenance  applications.  This  heuristic  is  designed  to 
have  the  following  quality  attributes:  (i)  ability  to  produce  an  approximate  solution  in  good  run 
time,  (ii)  ability  to  produce  approximate  solutions  that  are  close  to  exact  solutions,  (iii)  ability 
to  utilize  initial  solution  guesses  given  by  the  user  or  obtained  from  earlier  iterations,  (iv)  ability 
to  profitably  incorporate  new  results  on  standard  set-covering  and  related  combinatorial  problems. 

This  heuristic  consists  of  the  following  main  components:  (i)  a  discrete  line  search  (a  line 
search  wherein  the  function  is  evaluated  at  integral  points  only)  that  selects  combinations  of 
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sensor  mechanisms,  (ii)  a  random  selection  of  deployments  for  combinations  of  sensor 
mechanisms,  and  (iii)  a  "greedy"  procedure  for  computing  cover  value  for  any  given  deployment 
of  sensor  mechanisms. 

The  discrete  line  search  does  "skip  and  search"  operations  over  the  set  of  all  possible 
combinations  of  sensor  mechanisms.  First,  it  goes  through  all  k-subsets  (of  sensor  mechanisms), 
then  through  all  (k+1)- subsets,  and  so  on.  This  search  procedure  allows  the  user  to  input  a  guess 
of  how  many  sensor  mechanisms  may  suffice  to  cover  the  space  optimally,  that  is,  the  starting 
k,  but  this  is  not  required.  This  guess  is  used  to  initiate  the  search,  and  naturally  loses  its  effect 
if  it  is  misleading. 

If  a  selected  deployment  (of  a  given  combination  of  sensor  mechanisms)  does  not  cover 
the  space,  the  heuristic  makes  another  selection  of  deployment.  This  way,  each  combination  of 
sensor  mechanism  gets  a  number  of  "chances";  the  number  of  chances  given  to  each  combination 
is  determined  at  heuristic  initialization. 

Computational  experiments  indicate  that  the  "greedy  method"  for  computing  cover  value 
(cost  or  time)  is  a  good  choice.  One  ceui  show  that  a  related  divide-and-conquer  method  suited 
to  a  parallel  (processing)  computer  will  perform  even  better  than  that  greedy  method. 

Appendix  II  gives  other  details  of  this  heuristic  for  problem  (*).  It  also  includes  the 
source  code  of  a  computer  program  (in  C)  that  implements  the  heuristic.  This  program  is 
included  as  a  subroutine  in  the  larger  computer  program  used  for  the  computational  experiments 
discussed  in  section  3  of  this  report. 
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SECTION  3 


COMPUTATIONAL  EXPERIMENTS 

The  computational  experiments  reported  here  serve  two  purposes.  First,  they  are  designed 
to  demonstrate  that  our  procedure  for  problem  (*)  generates  solutions  close  to  optimal  solutions. 
Secondly,  they  are  used  to  examine  some  of  the  design  choices  that  define  the  procedure.  Such 
a  numerical  evaluation  of  design  choices  may  be  necessary  even  if  performance  guarantees  have 
already  been  given  analytically. 

3.1  Design  Choices 

The  main  design  choices  in  our  procedure  for  problem  (*)  are: 

i.  the  use  of  binary  code  to  "line-up"  all  combinations  of  sensor  mechanisms, 

ii.  the  restart  feature,  instead  of  a  one-pass  search, 

iii.  the  use  of  random  deployment,  for  any  given  combination  of  sensor  mechanisms, 
instead  of  an  orderly  enumeration  of  possibilities,  and 

iv.  the  use  of  a  greedy  method  for  computing  the  objective  function  values  (once  the 
deployment  is  determined). 

Alternatives  to  the  binary  code,  such  as  the  gray  code,  were  considered  before  the  binary 
code  was  chosen.  Computational  experiments  (using  a  program  written  in  BASIC  with  function 
plotting,  etc.)  put  the  binary  code  ahead  of  other  options  for  two  reasons.  First,  the  binary  code 
is  easy  to  explain.  Secondly,  the  function  plots  obtained  for  binary  codes  were  more  amenable 
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to  function  minimization  than  they  are  for  other  codes.  To  explain  this  apparent  superiority  of 
the  binary  code,  we  focused  on  k-terms  (that  is,  terms  of  the  binary  code  sequence  that  have 
exactly  k  I’s  in  their  digits),  and  discovered  some  orderliness  that  led  us  to  the  "restart"  idea. 

To  see  this  k-term  property  of  binary  code  sequences,  consider  the  k-terms  of  the 
following  sequence. 

0000,  0001,  0010,  0011,  0100,  0101,  0110,  0111, 

1000,  1001,  1010,  1011,  1100,  1101,  1110,  nil. 

The  sequences  of  k-terms  are 

0001,  0010,  0100,  1000, 

0011,  0101,  0110,  1001,  1010,  1100, 

0111,  1011,  1101,  1110, 

nil. 

These  suggested  to  us  that  it  might  be  advantageous  to  concentrate  the  line  search  on  one 
sequence  of  k-terms  at  a  time.  Hence,  the  idea  of  restart.  Computational  experiments  with  test 
problems  (described  below)  suggest  that  the  "restart"  feature  is  indeed  a  good  idea,  besides 
allowing  the  user  to  input  a  guess  about  how  many  sensor  mechanisms  may  suffice  to  cover  the 
space. 

With  regard  to  randomness  in  the  deployment  choice,  the  only  alternative  is  to  enumerate 
possibilities  on  the  basis  of  information  obtained  from  some  pre-processing  of  data  concerning 
deployment  options.  We  did  not  perform  computations  with  this  alternative  because  it  seems  too 
much  work  for  too  little  reward.  The  deployment  choice  aspect  has  significant  scope  for  future 
work.  However,  the  randomness  seems  to  have  no  adverse  effect  on  the  results  when  relatively 
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few  of  the  sensor  mechanisms  have  multiple  deployment  options. 

The  use  of  a  greedy  method  for  computing  objective  function  values  is  convenient  on 
account  of  speed.  Any  sensible  heuristic  for  the  standard  set-covering  problem  may  do  equally 
well.  A  "divide-and-conquer"  alternative  was  considered,  but  this  is  more  suited  to  a  parallel 
(processing)  computer. 

3.2  Hardware  and  Programming 

To  some  extent,  all  computational  experiments  reflect  hardware  and  programming.  Our 
computer  program  for  implementing  the  procedure  is  written  in  C  language  so  as  to  take 
advantage  of  C’s  efficacy  in  handling  complicated  data  structures.  We  had  presumed  that 
subdivisions  of  the  set  S  might  require  such  data  structures. 

Development  started  with  an  IBM  RT  PC  in  a  UNIX  environment,  but  we  switched  to 
an  IBM  AT  PC  (286)  and  an  IBM  PS/2  Model  80,  so  as  to  take  advantage  of  existing  utility 
packages.  The  source  code  included  here  as  an  Appendix  is  in  Microsoft  C,  using  function 
prototyping  and  other  features  of  ANSI  C.  Our  test  problems  are  of  the  size  that  will  run  with 
DOS.  It  is  a  simple  matter  to  make  the  little  adjustments  (function  declaration,  etc.)  that  will 
make  this  code  run  in  the  environment  of  other  operating  systems. 

Our  program  accepts  inputs  in  form  of  files  stored  in  "binary"  (for  space  use  efficiency). 
The  user  is  required  to  supply  the  file  name  when  the  file  is  to  be  read  or  written  into.  The  user 
also  selects  the  objective  function  -  cost  or  time.  Next,  the  user  is  asked  to  guess  (if  possible) 
how  many  sensor  mechanisms  may  suffice  to  cover  the  space. 
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3.3  Test  Problems 


Our  test  problems  are  not  of  the  most  general  kind.  They  have  been  designed  with  the 
surveillance  application  in  mind.  We  believe  that  the  surveillance  application  will  ordinarily 
have  (i)  many  pieces  in  the  cubical  subdivision  of  S  (Appendix  I),  (ii)  few  sensor  mechanisms, 
(iii)  most  sensor  mechanisms  having  exactly  one  deployment  option.  Test  problems  reflecting 
these  observations,  and  computed  solutions  are  given  in  Appendix  IV. 

Each  problem  has  six  sensor  mechanisms,  and  each  sensor  mechanism  has  at  most  two 
deployment  options.  The  coverage  span  of  each  sensor  mechanism  is  assumed  to  be  a  sphere, 
and  the  set  S  is  assumed  to  be  a  rectangular  volume.  In  our  experience,  when  computed  "cost" 
solution  is  not  optimal,  the  computed  "time"  solution  usually  gives  help.  In  any  case,  supplying 
"guesses"  always  helps. 

It  is  clear  from  an  inspection  of  those  test  problems  that  there  is  need  for  more 
computational  experiments,  especially  ones  involving  large-scale  problems,  the  kind  that  may 
require  a  computing  environment  such  as  that  provided  by  OS/2  or  UNIX. 
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DIRECTIONS  FOR  FUTURE  WORK 

Directions  for  further  research  related  to  problem  (*)  are  suggested  in  this  section.  Three 
classes  of  extensions  of  problem  (*)  will  be  described.  The  first  class  has  to  do  with  making 
problem  (*)  more  realistic  in  terms  of  resource  constraints,  for  example,  budget  constraint  and 
time  constraint.  The  second  class  is  about  how  existing  commitments  may  constrain  the  choice 
of  sensor  mechanisms  and  their  deployment  options.  The  third  class  of  extensions  is  concerned 
with  multiple-layer  surveillance  coverage  for  selected  portions  of  the  space  under  surveillance. 
In  addition,  further  computational  experiments  on  our  procedure  for  problem  (*)  will  continue 
to  lie  of  interest. 

4.1  Resource-Constrained  Extensions 

Each  time  an  instance  of  problem  (*)  is  solved,  the  objective  function  is  either  "cost"  or 
"time".  But,  in  the  real  world,  one  may  be  interested  in  both  time  and  cost  simultaneously.  It 
may  be  desirable  either  to  minimize  surveillance  rescheduling  cost,  with  a  limit  on  the  time  it 
takes  to  complete  the  rescheduling,  or  to  minimize  the  time  it  takes  to  complete  surveillance 
rescheduling,  with  a  limit  on  the  cost  of  doing  the  rescheduling. 

The  following  is  a  formal  statement  of  these  two  resource-constrained  extensions  of 
problem  (*): 

(1)  Problem  (*),  with  objective  function  replaced  by  minimize  subject  to  the 
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constraint  that 


Sjj  is  in  the  collection,  and  max{tjj)^T 
(2)  Problem  (*),  with  objective  function  replaced  by  minimize  max{tij}  subject  to  the 
constraint  that  S^j  is  in  the  collection  and  S  Cy  <.C 

Here,  T  is  the  time  limit,  and  C  is  the  cost  limit.  These  resource-constrained  problems 
are  "mixed",  because  each  problem  is  a  mixture  of  combinatorial  and  continuous  variable 
optimization  problems. 

4.2  Commitment-Constrained  Extensions 

In  the  real  world,  existing  contracts,  concerning  orders  and  supplies,  will  ordinarily 
constrain  new  decisions.  Therefore,  decisions  involved  in  rescheduling  surveillance  coverage 
may  be  subject  to  commitment  constraints.  Existing  commitments  on  orders  and  supplies  for 
surveillance  maintenance  may  be  reflected  in  problem  (*)  in  form  of  pre-determined  selection  of 
some  sensor  mechanisms  or  some  deployment  options.  Indeed,  the  standard  set-covering  problem 
may  be  regarded  as  an  extreme-case  member  of  this  class  of  extensions  of  problem  (*);  it  is  the 
case  where  each  surveillance  sensor  mechanism  is  allowed  exactly  one  deployment  option,  an 
option  probably  fixed  by  some  legal  contract. 

4.3  Multiple-Layer  Surveillance  Coverage 

It  is  possible  that  some  portions  of  the  space  under  surveillance  in  problem  (*)  are  so 
special  that  they  require  a  multiple-layer  surveillance  coverage.  For  example,  suppose  it  is 
necessary  to  ensure  that  portions  of  the  boundaries  of  the  space  be  covered  by  at  least  2  sensor 
mechanisms.  This  may  be  required  to  satisfy  some  reliability  requirements.  We  would  refer  to 
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this  as  a  2-layer  coverage  problem. 

This  gives  rise  to  a  k-layer  extension  of  problem  (*).  Problem  (*)  is  obviously  a  1 -layer 
coverage  problem.  Now,  if  we  allow  portions  of  the  space  to  have  0-layer  coverage,  and  other 
portions  to  have  positive-layer  coverage,  then  we  have  a  very  practical  extension  of  problem  (*). 
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SUBDIVIDING  3-MANIFOLDS  FOR  SURVEILLANCE  MAINTENANCE 

by 

S.  Awoniyi  &  L.  Frair 

Abstract 

This  report  contains  a  procedure  for  checking  whether  a  given  set  of  surveillance  sensor 
mechanisms  is  sufficient  to  cover  a  given  airspace,  say  S,  under  surveillance.  This  procedure 
consists  of  using  information  about  the  size  of  the  smallest  potential  intruder  (into  S)  to  subdivide 
S  into  pieces  that  fit  together  suitably.  The  vertices  of  these  pieces  constitute  a  finite  set  of 
points  that  the  sensor  mechanisms  must  cover.  Hence,  this  procedure  essentially  replaces  S  with 
a  finite  set;  then,  checking  whether  S  is  covered  by  the  sensor  mechanisms  becomes  a 
straightforward  task  on  a  digital  computer.  Also  included  (as  Appendix  II)  in  this  report  is  a 
menu-driven  computer  program  (in  C  programming  language)  that  implements  that  procedure; 
it  subdivides  S  and  denumerates  (assigns  numbers  to)  its  pieces,  in  a  manner  that  helps  record¬ 
keeping.  This  computer  program  will  be  incorporated  into  the  computer  codes  that  will  be 
delivered  to  the  US  Army  Strategic  Defense  Command  under  the  research  contract  DASG-60-90- 
C-0142.  In  the  mean  time,  this  computer  program  may  be  used  for  experimentation  and 
demonstrations. 
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1.  Objective  of  Results 

The  results  given  in  this  report  address  aspects  of  the  generalized  set-covering  problem 
described  in  the  research  grant  DASG-60-90-C-0142,  sponsored  by  the  United  States  Army 
Strategic  Defense  Command.  The  problem  is  recalled  here: 

(*)Giyen:  S,  a  subset  of  the  3-dimensional  space,  R^;  subsets  Sjj,  i=L..,  m,  j=l,...,n, 
that  cover  S  in  the  sense  that  S  is  the  union  of  all  the  S^’s;  two  numbers  t^j  and 
Cy  that  will  be  interpreted  as  "time"  and  "cost"  respectively.  Required:  To 
determine  an  optimal  collection  of  the  S^’s  that  cover  S,  with  no  more  than  one 
Sij  for  each  i,  and  optimality  defined  in  terms  of  minimizing  max  {tjj:Sij  is  in  the 
collection}  or  minimizing  the  total  costs  of  the  Sy’s  in  the  collection. 

This  set-covering  problem  shall  henceforth  be  referred  to  as  problem  (*). 

Any  method  for  solving  problem  (*)  must  include  a  procedure  for  checking  whether  the 
set  S  is  equal  to  the  union  of  some  of  the  Sj^’s.  If  the  set  S  were  a  finite  set,  then  it  would  be 
a  straightforward  matter  to  do  this  checking  on  a  digital  computer. 

In  the  surveillance  applications  of  problem  (*),  the  set  S  is  not  finite;  S  is  a  3-dimensional 
manifold  (that  is,  a  full-dimensional  volume  in  R^).  S  represents  the  space  under  surveillance, 
whereas  Sy  represents  the  portion  of  S  that  may  be  covered  by  the  i-th  surveillance  sensor 
mechanism  functioning  in  its  j-th  optional  positioning. 

If  S  is  a  3-manifold  without  any  special  features,  to  determine  whether  S  is  equal  to  the 
union  of  some  of  the  Sy’s  is  a  very  hard  task  on  a  digital  computer.  But  the  surveillance 
maintenance  situation  allows  us  to  regard  S  as  "variably  finite"  (this  term  will  be  fully  explained 
in  Section  2).  Utilizing  this  special  feature  of  S,  we  give  in  this  report  a  reasonable  procedure 
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for  comparing  S  and  the  union  of  any  collection  of  the  Sy’s. 

This  procedure  consists  of  (a)  subdividing  S  into  a  finite  number  of  pieces  of  suitable 
sizes,  with  each  piece  "contained"  in  at  least  one  of  the  Sy’s,  and  (b)  denumerating  (numbering) 
the  pieces  in  a  way  that  will  facilitate  efficient  record-4ceeping.  These  results  essentially  replace 
S  with  a  finite  set  { <r,  ^ .  cTj, } ,  where*3j^  is  the  k-th  piece  of  the  subdivision,  and  p  the  total 
number  of  such  pieces.  For  overall  efficiency,  it  is  desirable  to  make  p  as  small  as  possible, 
while,  at  the  same  time,  ensuring  that  it  does  not  take  much  effort  to  check  whether  a  piece  is 
covered  by  an  Sy.  Section  2  of  this  report  gives  details  of  these  results.  A  computer  code  (in 
C  language)  implementing  this  procedure  is  included  as  Appendix  II. 

2.  Description  of  Results 

This  section  describes  our  procedure  for  comparing  S  and  unions  of  the  S^’s.  We  shall 
begin  by  making  some  problem-reducing  observations  on  the  surveillance  maintenance  situation. 
These  observations  are  instrumental  to  our  replacing  S  with  a  finite  set,  thereby  making  the  set 
comparison  task  a  straightforward  one. 

(*)Observations:  (i)  S  is  compact,  since  surveillance  over  an  unbounded  space  is  not 
practicable;  (ii)  surveillance  sensor  mechanisms  do  not  have  to  cover  every  point  of  S; 
they  only  need  to  cover  points  in  S  that  are  not  too  far  apart  relative  to  the  size  of  the 
smallest  potential  intruder  into  S. 

To  substantiate  the  claim  contained  in  the  second  part  of  those  observations,  let  us 
suppose  S  is  a  rectangular  volume,  and  the  smallest  potential  intruder  into  S  is  of  size  « .  Let 
S  be  subdivided  into  polyhedral  pieces  (that  is,  pieces  with  straight  boundaries)  {a',,...,crp),  each 
with  diameter  less  than  o(  (the  diameter  of  a  set  is  the  length  of  the  longest  straightline 
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contained  in  the  set).  The  compactness  of  S  guarantees  the  existence  of  such  polyhedra  pieces. 
Obviously,  the  surveillance  sensor  mechanisms  only  need  to  cover  the  vertices  of  thee"’s,  since 
that  guarantees  that  any  intruder  will  be  detected.  Figure  1  gives  an  illustration  in  of  possible 
subdivisions  of  manifolds.  Thus,  S  is  effectively  replaced  with  a  finite  set,  the  set  of  all  vertices 
of  the  pieces  of  a  subdivision  of  S.  Hence,  by  virtue  of  the  observations  (*),  S  may  be  regarded 
as  a  "variably  finite  set".  The  qualifier  "variably"  is  necessary  because  there  are  many  different 
wrys  of  subdividing  S  into  pieces  giving  rise  to  many  different  finite  sets  for  replacing 

S. 

As  one  can  see  from  Figure  1,  some  subdivisions  of  S  are  better  than  others.  In  choosing 
a  subdivision  of  S,  it  is  desirable  to 

(i)  minimize  the  number  of  (subdivision)  pieces;  Figures  1(a)  and  1(c)  show  examples 
of  subdivisions  with  relatively  too  many  pieces; 

(ii)  avoid  repeating  a  vertex  treatment  in  assigning  (subdivision)  pieces  to  the  Sjj’s;  such 
a  repetition  will  occur  if  there  are  more  pieces  than  vertices;  Figure  1(a)  illustrates 
this  point; 

(iii)  ensure  effective  denumeration  of  pieces  (effective  in  the  sense  that  each  piece  number 

automatically  identifies  the  corresponding  piece,  and  vice  versa),  so  as  to  facilitate 
efficient  record  keeping;  in  this  regard,  the  subdivision  of  Figure  1(b)  is  superior  to 
those  of  Figures  1(a)  and  1(c). 

Evidently,  the  usual  performance  criteria  for  subdivisions  in  homotopy  algorithms  [5]  , 
such  as  "directional  density",  do  not  apply  here  because  subdivision  traversing  is  not  of  interest 
in  the  context  of  surveillance  maintenance.  Appendix  I  gives  a  formalization  of  criteria  (i)  and 
(ii)  above. 
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On  the  basis  of  the  three  criteria  listed  above,  the  subdivision  of  Figure  1(b)  is  clearly 
superior  to  those  of  Figures  1(a)  and  1(c).  In  fact,  in  R^,  when  one  considers  using  k-sided 
pieces,  for  k=3,4,...,  it  is  not  hard  to  see  that  the  subdivision  in  Figure  1(b)  is  superior  to  all 
others.  In  R^,  analogues  of  the  subdivision  of  Figure  1(c)  fail  our  criteria  because  they  involve 
too  many  vertices  and  are  too  complicated  for  effective  denumeration;  such  subdivisions  are 
never  used  in  the  literature  of  homotopy  algorithms  and  finite  element  methods  [2]. 

Therefore,  the  feasible  options  reduce  to  choosing  between  3-dimensional  analogues  of 
the  subdivision  of  Figure  1(a),  which  are  known  as  "triangulations",  and  3-dimensional  analogues 
of  the  subdivision  of  figure  1(b),  which  we  shall  term  "cubical  subdivisions".  Cubical 
subdivisions  cover  R’  with  identical  cubes  or  rectangular  blocks,  whereas  triangulations  cover  R’ 
with  tetrahedra.  Now,  triangulations  fail  our  criteria  on  account  of  having  more  pieces  (called 
simplices)  than  vertices  (so  vertex  treatment  repetition  will  occur),  and  because  they  are  difficult 
to  denumerate  effectively.  In  contrast,  the  cubical  subdivision  has  relatively  few  pieces  (always 
less  than  the  number  of  vertices),  and  is  amenable  to  easy,  effective  denumeration.  Indeed,  for 
the  cubical  subdivision^  only  one  vertex  of  each  piece  (its  origin)  needs  to  be  assigned  to  the  S^j’s, 
because  the  other  vertices  (of  the  piece)  are  origins  for  other  pieces,  unless  the  piece  happens  to 
lie  on  some  of  the  boundary  of  S.  Figure  2  illustrates  this  point  about  the  cubical  subdivision. 
Hence,  we  have  a  case  for  choosing  the  cubical  subdivision  for  our  purposes  here. 

But  to  conclude  the  case  for  the  cubical  subdivision,  it  is  necessary  to  consider  the  amount 
of  record-keeping  that  must  be  done  at  the  boundary  of  S,  taking  into  account  the  types  of  S  that 
may  arise  in  the  surveillance  maintenance  situation  (consideration  of  the  interior  of  S  having  been 
subsumed  in  the  foregoing  comparison  of  subdivisions  of  S).  Now,  since,  by  definition,  S 
represents  an  airspace,  one  may  assume  that  S  is  essentially  a  3-manifold  of  the  type  shown  in 
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S  =  3  X  4  squarej  c<= 
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Figure  3(a).  As  the  manifold  of  Figure  3(a)  is  clearly  homeomorphic  to  that  in  Figure  3(b),  in 
a  topological  sense,  one  may  then  assume  that  S  is  essentially  a  rectangular  volume,  without  any 
loss  of  generality.  Hence,  S  indeed  possesses  the  features  that  have  been  demonstrated  ( in  the 
foregoing)  as  reinforcing  the  cubical  subdivision’s  superiority.  This  concludes  our  case  for 
choosing  the  cubical  subdivision. 

The  foregoing  analysis  is  the  basis  for  a  computer  program  that  is  included  here  as  an 
appendix.  This  computer  program  gives  a  cubical  subdivision  of  S,  and  denumerates  its  pieces. 
It  will  be  included  as  a  component  of  the  final  computer  implementation  of  our  solution  of 
problem  (*).  Presently,  this  computer  program  may  be  used  in  computational  experiments  and 
demonstrations  on  the  cubical  subdivision. 

The  computer  program  is  menu-driven  and  user-friendly.  First,  it  requires  the  user  to 
supply  the  dimensions  (length,  breadth,  height)  of  the  rectangular  volume  S.  Then  the  user 
makes  a  selection  from  the  following  menu: 

A  -  To  Compute  Number  Assigned  to  a  Piece  ^ 

B  -  To  Compute  the  Piece  for  a  Given  Number 

C  -  To  Return  To  Operation  System 

If  A  is  selected,  then  the  user  is  asked  to  supply  the  coordinates  of  the  piece’s  origin;  the 
computer  program  then  returns  the  piece’s  number.  If  B  is  selected,  then  the  user  is  asked  to  give 
the  number,  the  program  then  returns  the  origin  of  the  corresponding  piece,  unless  the  given 
number  is  too  big  to  be  a  piece  number.  In  the  case  that  the  number  is  too  big,  the  program 
returns  an  error  message. 


8 


3  (a)  S  =  Airspace 


3  (b)  S  After  Honeonorphic 
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Figure  3i  The  Form  of  S 
In  Problem  (*) 
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3.  Concluding  Remarks 

Besides  allowing  us  to  compare  S  to  collections  of  the  Sy’s,  replacing  S  with  a  finite  set 
also  transforms  problem  (*)  into  a  generalized  discrete  set-covering  problem  that  is  amenable  to 
a  host  of  discrete  optimization  strategies  and  heuristics.  These  strategies  include  divide-and- 
conquer,  discrete  line  search,  decomposition,  and  enumeration. 

To  employ  these  strategies  to  solve  problem  (*),  we  shall  use  the  computer  program 
included  here  as  a  "function"  (or  subroutine).  This  function  will  be  called  for  doing  "initial  set¬ 
up"  (implicitly  replacing  S  and  the  Sjj’s  with  finite  sets),  and  for  checking  for  feasibility  when 
collections  of  the  Sy’s  are  examined. 

Since  each  Sy  will  contain  relatively  few  of  the  pieces  (of  the  cubical  subdivision  of  S), 
Sy  may  be  represented  as  a  "linked  list".  The  C  programming  language  is  ideal  for  such  a  data 
structure. 

Now,  our  analysis  (in  Section  2)  that  resulted  in  the  choice  of  the  cubical  subdivision  is 
relatively  informal.  A  more  formal  analysis  resulting  in  the  same  conclusion  is  included  here  as 
Appendix  I. 
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APPENDIX  I 

A  FORMAL  CASE  FOR  THE  CUBICAL  SUBDIVISION 
Here,  we  present  in  a  more  rigorous  version  of  our  case  for  the  cubical  subdivision. 
We  shall  also  indicate  how  this  analysis  may  be  extended  to  R’.  First,  a  number  of  concepts  will 
be  defined. 

Deflnition:  For  any  polygon,  <r,  in  R^,  the  length  of  the  longest  side  shall  be  called  the 
polygon’s  surveillance  diameter,  and  will  be  denoted  by  sdiam(ff). 

(T 

The  number  sdiam(<i)  is  the  upper  bound  on  the  size  of  intruders  that  may  enter^undetected  by 
surveillance  sensors  that  cover  the  vertices  of  tr  only.  Our  next  definition  is  a  performance 
measure  (for  subdivisions)  summarizing  two  of  the  three  criteria  given  in  Section  2. 
Deflnition:  For  any  n-sided  polygon,  o',  in  R^  with  area  A,  the  number  A/n  shall  be  called  the 
polygon’s  index  of  surveillance  subdivision  efficacy,  and  will  be  denoted  by  isse(o*). 

I^for  a  polygon  ajsdiam(tr)  is  acceptably  small  and  isseCo")  is  relatively  large,  then 
cr  is  a  good  candidate  for  use  as  a  subdivision’s  basic  building  block. 

The  following  lemma  follows  easily  from  the  definitions  above: 

Lemma:  (i)  If  o'  is  an  equilateral  triangle  with  side  C{ ,  then 
sdiamCc^  =  and  isseCcr)  =  <x74'5’ 

(ii)  If  cr  is  a  square  of  side  length  a ,  then 
sdiam(e)  =  o<  and  isse((j)  =  o^V4 

Hence,  for  surveillance  purposes  in  R^,  it  is  better  to  subdivide  S  into  squares  rather  than 
equilateral  triangles,  especially  if  S  happens  to  be  a  rectangular  area.  However,  that  does  not  yet 
indicate  that  squares  are  preferable  to  all  triangles;  this  is  the  object  of  the  following  proposition. 
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Proposition:  For  any  triangle,  cr,  in  with  sdiam(cT)  =  oC  , 

isse(cr)  < 

Proof  (outline):  It  suffices  to  show  that  die  area  of  the  triangle  o' does  not  exceed 73ocV4  (the  area 

of  coaesponding  equilateral  triangle).  Since  sdiam(<r)  =  the  longest  side  ofiris  of  length 

With  a  simple  construction,  one  can  enlarge  <r  into  an  isosceles  triangle  with  two  sides  equal  -to  k  . 

Then,  using  an  elementary  optimization  technique,  one  can  show  that  if  CT  is  isosceles  with 
sdiam(<^  =  then  the  area  of  o'  does  not  exceed  that  of  an  equilateral  triangle  with  side  length 
equal  ta  pc . 

Thus,  we  have  formally  demonstrated  that  the  cubical  subdivision  of  S  is  superior  to 
triangulations  of  S  in  R^,  even  without  using  the  fact  that  S  is  a  rectangular  area.  This  conclusion 
also  holds  true  in  R’,  but  the  analysis  is  a  lot  longer.  In  R^,  the  tetrahedron  replaces  the  triangle 
of  R^  and  the  cube  replaces  the  square.  Of  course,  the  definitions  of  sdiam(T')  and  isse  C'J') 
must  be  recast  *-  "facet"  replacing  "line",  "volume"  replacing  "area",  etc. 

Obviously,  we  have  not  considered  k-sided  polygons  in  R^,  for  k=5,6,...  This  is  because 
such  polygons  are  too  complicated  (when  interpreted  for  R^)  to  allow  sensible  denumeration  of 
subdivision  pieces.  They  also  fit  poorly  on  the  boundary  of  S,  especially  when  S  is  a  rectangular 
volume. 
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1  of  2 

*  CUBICAL  SUBDIVISION  OF  3-MANIFOLDS  FOR  SURVEILLANCE  MAINTENANCE:  */ 

*  WE  ASSUME,  WITHOUT  LOSS  OF  GENERALITY,  THAT  THE  MANIFOLD  IS  A  */ 

*  RECTANGULAR  3-DIMENSIONAL  VOLUME;  THE  MANIFOLD  IS  SUBDIVIDED  INTO*/ 

*  UNIT  CUBES,  AND  THE  CUBES  ARE  DENUMERATED  FOR  EFFICIENT  RECORD  */ 

*  KEEPING;  THE  USER  SHALL  SUPPLY  THE  MANIFOLD'S  DIMENSIONS  (LENGTH,*/ 

*  BREADTH,  AND  HEIGHT),  AND  THEREAFTER  SELECT  FROM  A  MENU  THAT  */ 

*  GIVES  A  CUBE'S  NUMBER  WHEN  THE  CUBE'S  ORIGIN  COORDINATES  ARE  */ 

*  SUPPLIED,  AND  GIVES  A  CUBE'S  ORIGIN  COORDINATES  WHEN  THE  CUBE'S  */ 

*  NUMBER  IS  SUPPLIED.  */ 

oi d  pi snum{ voi d ) ; 
old  nunipls(void); 
nt  11,  bb ,  hh ; 
ain{ ) 

{ 

char  ss; 

printf("  THIS  PROGRAM  GIVES  A  CUBICAL  SUBDIVISION  OF  A\n  "); 
printf("  3-DIMENSIONAL  RECTANGULAR  BLOCK\n\n  "); 

pri nt f ( "Pi  ease  ENTER  the  dimensions  of  the  block,  ptqsr  -->  "); 
scanf { "%d: %d: %d" ,  &11,  &bb,  &hh); 
start : 

pri  ntf { "\n\n" )  ; 

printf ( "\n  ENTER  \n"  )  ; 

printf{"\n  A  --  To  Compute  Number  Assigned  to  a  Piece\n"); 

printf ("\n  B  --  To  Compute  the  Piece  for  a  Given  Number\n"); 

printf("\n  C  --  To  Rr  ti  -'n  to  Dds\n"); 

printf("\n  YOUR  SELECTION  r-->  "); 

ss  =  getcheO;  pri  ntf  (  "\n\n"  )  ; 

switch(ss) 

{ 

case  ' A '  : 

pi snum ( ) ;  break ; 
case  '  a  ' : 

pi snum { )  ;  break  ; 
case  ' B '  : 

numpi s ( ) ;  break ; 
case  ' b  '  : 

numpi s ( ) ;  break  ; 
case  'C'l 
exi t (0 )  ; 
c  a  s  e  '  c  '  : 

exit(0); 
def aul t ; 

goto  start; 

} 

pri nt  f ( "\n\n " ) ; 

printf  (  " - \n\n\n"  )  ; 

pri ntf ( "Press  any  key  to  continue"); 
s s  »  getch ( ) ; 
goto  start; 

} 
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2  of  2 


/*  THIS  PROCEDURE  COMPUTES  THE  NUMBER  ASSIGNED  TO  ANY  */ 

/*  PIECE  OF  THE  CUBIC  SUBDIVISION  OF  THE  BLOCK,  /nu  ARE  */ 

/*  REQUIRED  TO  INPUT  THE  ORIGIN  OF  THE  PIECE. \n\n»/ 

void  pi snum ( voi d ) 

Int  XX,  yy,  zz,  temp; 

start  1 :  _  _ 

pri nt f ( "Pi  ease  ENTER  the  origin  coordinates,  x:y:z  -->  "); 
scanf ( "*ds %dj %d"  ,  &xx,  &yy,  &2z); 
if(xx>ll  II  yy>bb  ||  2z>hh) 

{printf  ("\n\n!''!''!'‘!''!^!''!''!''INPUT  ERROR!  TRY  AGAIN  PLEASE\n\n\n "  )  ; 
goto  startl;  } 

temp  =  ( 1 1 *bb*zz ) +( 1 1 *yy ) +XX+1 ; 

pri ntf ( "\nFor  the  piece  with  origin  {%d,%d,%d),  ",  xx,  yy,  zz); 
printf("the  corresponding  number  is  %4d  \n",  temp); 


/*  THIS  PROCEDURE  COMPUTES  THE  ORIGIN  OF  THE  PIECE  */ 

/*  CORRESPONDING  TO  ANY  GIVEN  NUMBER*/ 

voi d  numpi s (voi d ) 

int  num,  xx,  yy,  zz,  temp;  >_ 

St art 2: 

pri ntf ( "PI  ease  ENTER  the  number  -->  "); 
scanf("%d",  &num);  pri nt f ( "\n " ) ; 
if  ( num>( 1 1 *bb*hh )  ) 

{printf  ("\n\n!^!^l''!''!^!''l''!''INPUT  ERROR!  TRY  AGAIN  PLEASE\n\n  "  )  ; 
goto  start2; 

} 

num--l;  temp  =  ll*bb; 

zz  »  (  num  -  (numttemp)  )/temp;  num  zz*temp; 
yy  »  {  num  -  (numtll)  )/ll; 

XX  ■  num  -  yy* 1 1  ; 

printf("The  corresponding  piece  has  {%d,td,%d)  ",  xx,  yy,  zz); 

printf{"as  origin  coordi nates\n\n" ) ; 
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A  RESTART  PROCEDURE  FOR  A  GENERALIZED 
SET-COVERING  PROBLEM  IN  SURVEILLANCE  MAINTENANCE 


by 

Samuel  Awroniyl  &  Lesber  Frair 


Abstract :  We  describe  a  heuristic  procedure  for  the  generalized 
set-covering  problem  defined  in  our  research  contract  DASG 
60-90-C-0142  with  the  United  States  Army  Strategic  Defense  Command, 
Huntsville,  Alabama.  This  heuristic  is  not  a  variant  of  any 
existing  procedure,  as  the  problem  itself  has  not  appeared  in  the 
literature.  The  quality  attributes  of  this  heuristic  include 

(i)  ability  to  produce  an  approximate  solution  in  good  run  time, 

(ii)  ability  to  produce  approximate  solutions  that  are  close  to 
exact  solutions,  (iii)  ability  to  improve  upon  initial  solution 
guesses  given  by  the  user  or  obtained  from  earlier  iterations, 
(iv)  ability  to  profitably  incorporate  new  results  on  set-covering 
and  related  combinatorial  problems.  The  main  components  of  this 
heuristic  are  a  discrete  line  search  with  restart,  a  random 
selection  operation,  and  a  greedy  procedure  on  a  subproblem.  A 
computer  program  (in  C  language)  implementing  this  heuristic  is 
included  as  an  Appendix.  The  next  report  on  this  research  contract 
(which  is  also  the  final  report)  will  give  a  full  account  of  our 
computational  experiments  with  this  heuristic. 
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1 .  Introduction 

This  is  the  second  of  three  technical  reports  fulfilling  a 
requirement  of  our  research  contract  DASG  60-90-C-0142  with  the  US 
Army  Strategic  Defense  Command,  Huntsville,  Alabama.  The  goal  of 
this  research  contract  is  to  design  solution  procedures  for  the 
generalized  set-covering  problem  described  in  (*)  below. 

(*) Given:  S,  a  subset  of  the  3-dimensional  space,  R^;  subsets 
Sij/  i=l/  .  .  .,ni,  j=l,  .  .  .,n,  that  cover  S  in  the  sense  that  S 
is  the  union  of  all  the  Sj.j's;  two  numbers  t^^j  and  c^j  that 
will  be  interpreted  as  "time"  and  "cost"  respectively. 
Required:  To  determine  an  optimal  collection  of  the  S^^'s 
that  cover  S,  with  no  more  than  one  Si^  for  each  i, 
and  optimality  defined  in  terms  of  minimizing 
max{tij  :  S^^  is  in  the  collection}  or  minimizing  the  total 
cost  of  the  Sij's  in  the  collection. 

This  set-covering  problem  shall  be  referred  to  as  problem  (*)  . 
Problem  (*)  is  a  representation  for  a  class  of  surveillance 
maintenance  problems  including  the  problem  of  restoring 
surveillance  coverage  for  a  system  of  surveillance  sensor 
mechanisms  after  a  sensor  mechanism  (in  the  system)  fails.  The  set 
S  represents  the  space  under  surveillance;  each  subset 
represents  the  portion  of  S  covered  by  the  i-th  surveillance  sensor 
mechanism  functioning  in  its  j-th  optional  deployment/positioning; 
each  Cij  represents  the  cost  of  the  j-th  deployment  of  the  i-th 
sensor  mechanism;  t^^  represents  the  time  required  for  the  i-th 
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sensor  mechanism  to  be  placed  in  its  j-tb  optional  deployment. 
Selecting  a  minimum-cost  set  of  S^j's  covering  S  corresponds  to 
restoring  surveillance  coverage  at  minimum  cost,  whereas  selecting 
a  minimax-time  set  of  S^^'s  covering  S  corresponds  to  restoring 
surveillance  coverage  within  the  shortest  time  possible. 

Our  first  technical  report  is  entitled  "Subdividing  3- 
Manifolds  For  Surveillance  Maintenance",  and  our  third  technical 
report  shall  consist  of  a  summary  together  with  a  user-friendly 
computer  implementation  of  the  results  presented  in  the  first  two 
technical  reports .  The  third  report  will  also  give  indications 
of  further  work  that  should  be  done  in  order  to  profitably  adopt 
the  results  of  these  reports  in  the  real  world. 

This  second  technical  report  describes  a  practical  solution 
procedure  for  problem  (*)  .  To  be  able  to  see  the  origins  of 
certain  features  of  this  procedure,  it  is  necessary  to  understand 
the  combinatorial  nature  of  problem  (*)  .  Suppose  there  are  d^ 
deployment  possibilities  for  the  i-th  sensor  mechanism.  Then, 
there  are  (di+1)  (dj+l)  .  .  .  (dji+l)  valid  combinations  of  the  S^j's  that 
must  be  considered  explicitly  or  implicitly  by  any  procedure  that 
computes  an  exact  solution  to  problem  (*)  .  If  di=n,  each  i,  then 
the  niJmber  of  combinations  of  the  s  is  (n+l)“. 

Suppose,  also,  that  p  computational  steps  are  required  to 
compute  the  (objective  function)  value  of  each  of  those 
combinations  of  S^^'s.  Then,  the  maximum  number  of  computational 
steps  needed  to  compute  an  exact  solution  to  problem  (*)  is 
(di+1)  (dj+l)  .  .  .  (d.+l)p.  Since  we  are  assuming  that  it  does  not 
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depend  on  the  S^^'s,  the  number  p  ordinarily  depends  on  the  size  of 
the  space  S;  therefore,  we  shall  henceforth  write  it  as  p(|S|). 

We  shall  refer  to  the  number  (di+1)  (dj+l)  •  •  •  (d»+l)  p(|Sl)  as 
the  computational  upper  bound  (c.u.b.)  for  problem  (*) .  The  c.u.b. 
has  some  practical  implication.  To  explain  this  practical 
implication  of  the  c.u.b.,  let  us  consider  a  well-known  special 
case  of  problem  (*) ,  namely  the  standard  set-covering  problem. 
When  di=l,  all  i,  (that  is,  each  sensor  mechanism  has  exactly  one 
deployment  option) ,  problem  (*)  reduces  to  the  standard  set¬ 
covering  problem  (by  virtue  of  results  presented  in  our  first 
technical  report),  and  the  c.u.b.  becomes  2”p(|S|).  For  instance, 
for  p(|S|)  >  1  and  m  =  50,  the  c.u.b.  for  the  set-covering  problem 
is  at  least  2®**,  and  this  number  of  computational  steps  will  take 
a  fast  digital  computer  doing  one  operation  per  micro-second  about 
35  years  to  complete  (see  pp  141  of  [3]  for  more  on  computational 
speeds  and  their  computer  time  requirements) . 

One  clear  conclusion  from  the  foregoing  observations 
concerning  c.u.b.  is  that  if  a  procedure  must  compute  an  exact 
solution  for  an  instance  of  problem  (*) ,  then,  unless  m  is  suitably 
small,  that  procedure  must  keep  p(|S|)  small  enough  and  avoid 
explicit  consideration  of  most  of  the  combinations  of  the  S^j's. 
Now,  this  has  a  bearing  on  the  theory  of  algorithms .  It  is  well 
known  in  the  theory  of  algorithms  that  the  set-covering  problem  is 
NP-complete  [1]  [7];  this  means  that  it  is  highly  unlikely  that  one 
can  ever  find  a  practical  solution  procedure  (for  the  set-covering 
problem)  that  keeps  p(|S|)  small  enough  and,  at  the  same  time. 
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avoids  explicit  consideration  of  a  good  number  of  those 
combinations  of  the  S^^'s.  Our  third  technical  report  expatiates 
further  on  this. 

Eence,  in  solving  problem  (*)  in  the  real  world,  unless  S  and 
m  are  suitably  small,  approximate  solutions  should  be  the  goal; 
seeking  an  exact  solution  might  take  too  long.  Such  approximate 
solutions  should  not  be  arbitrary,  but  must  satisfy  sensible 
specifications  depending  on  the  application.  In  section  2  of  this 
report,  we  describe  attributes  that  we  consider  desirable  for  an 
approximate  solution  procedure  for  problem  (*),  in  view  of  the 
applications  described  above.  We  shall  follow  standard 
terminology,  and  refer  to  approximate  solution  procedures  as 
heuristic  procedures  or  heuristics. 

In  section  3,  our  heuristic  procedure  for  problem  (*)  is 
described.  Section  4  describes  our  computer  program  (in  C)  for 
implementing  this  heuristic,  and  the  program' s  source  code  is 
included  as  an  Appendix.  Our  next  technical  report  will  give  a 
complete  account  of  our  computational  experiments. 

2 .  Quality  Attributes  of  Heuristics  for  Problem  (*) 

Since  heuristics  are,  by  definition,  designed  to  produce 
approximate  solutions,  rather  than  exact  solutions,  it  is  necessary 
to  require  of  a  heuristic  certain  attributes  that  constitute  what 
may  be  termed  "good  quality".  (Of  course,  a  heuristic  may  very 
well  produce  exact  solutions  most  of  the  times)  .  For  any 
heuristic,  desirable  quality  attributes  would  ordinarily  be 
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determined  by  the  circumstances  of  particular  real  world 
applications.  Hence,  for  the  same  problem  model,  different 
heuristics  would  be  suited  to  different  real  world  applications . 

In  this  report,  we  shall  regard  problem  (*)  as  having  two  broad 
application  types,  namely,  the  surveillance  maintenance  application 
with  cost  minimization  as  objective,  and  the  surveillance 
maintenance  application  with  maximum-time  minimization  as 
objective.  It  turns  out  that  these  two  application  types  require 
slightly  different  heuristic  attributes. 

For  the  cost  minimization  version  of  problem  (*) ,  it  is 
desirable  that  a  heuristic  be  able  to  produce  approximate  solutions 
frequently  close  to  exact  solutions.  The  heuristic's  speed  is  of 
lesser  significance  than  the  quality  of  approximate  solutions 
generated.  However,  if  a  heuristic  has  a  way  of  improving  upon 
initial  solution  guesses,  then  good  speed  would  ordinarily 
translate  into  approximate  solution  of  good  quality  after  repeated 
iterations . 

For  the  maxim\am-time  minimization  version  of  problem  (*) ,  it 
may  be  desirable  that  a  heuristic  be  fast,  while  at  the  same  time 
producing  an  approximate  solution  of  good  quality.  (Recall  that 
the  maximum-time  minimization  version  is  a  representation  for  the 
problem  of  restoring  surveillance  cover  as  quickly  as  possible) . 
As  in  the  case  of  cost  minimization,  if  a  heuristic  can  produce  an 
approximate  solution  that  improves  upon  an  initial  solution  guess, 
then  good  speed  should  result  in  a  good  final  approximate  solution. 
Good  speed  together  with  quality  solutions  would  enhance  the 
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usefulness  of  flie  heuristic,  especially  in  emergency  situations . 

In  general,  a  heuristic  for  problem  (*) ,  irrespective  of  type 
of  objective  function,  should  be  designed  to  be  able  to  incorporate 
future  developments  on  set-covering  and  related  combinatorial 
optimization  problems .  This  is  an  especially  desirable  heuristic 
attribute  because  problem  (*)  is  practical,  NP-complete,  new  (in 
the  literature) ,  and  likely  to  generate  significant,  new  research 
activities  in  the  near  future. 

Finally,  as  the  question  of  surveillance  maintenance  concerns 
practical  system  security,  usually  involving  the  protection  of 
precious  human  lives  and  property,  it  is  desirable  that  a  heuristic 
for  problem  (*)  be  capable  of  being  described  in  terms  that  may  be 
easily  comprehended  by  persons  who  have  day-to-day  responsibility 
for  such  surveillance  systems .  Such  clarity  should  enhance  the 
heuristic's  chances  of  being  adopted  in  real-world  situations. 

3 .  Our  Heuristic  For  Problem  (*) 

We  describe  in  this  section  our  heuristic  for  problem  (*)  . 
This  heuristic  is  designed  to  have  the  attributes  described  in 
section  2.  In  summary,  the  attributes  are:  (i)  ability  to  produce 
an  approximate  solution  in  good  run  time  (that  is,  good 
computational  complexity) ,  (ii)  ability  to  produce  approximate 
solutions  that  are  close  enough  to  exact  solutions,  (iii)  ability 
to  improve  upon  initial  solution  guesses  given  by  the  user  or 
obtained  from  earlier  iterations,  (iv)  ability  to  profitably 
incorporate  new  results  on  set-covering  and  related 
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combinatorial  problems,  (v)  ease  of  comprehension,  so  as  to  secure 
user's  trust. 

Current  literature  contains  nothing  about  problem  (*)  in  its 
general  form.  Results  on  the  set -covering  problem  [4-]  [7]  are  all 
that  exist  in  current  literature  relating  directly  to  problem  (*) . 
Therefore,  our  heuristic  for  problem  (*)  is  not  a  variant  of  any 
existing  procedure.  The  following  is  an  informal,  application- 
oriented  description  of  the  heuristic. 

Heuristic 

Step  1:  Systematically  skip  and  search  through  (a  denxameration 
of)  subsets  of  the  given  set  of  sensor  mechanisms  so  as  to  find  a 
combination  of  sensor  mechanisms  possibly  "better”  than  the 
current  combination,  if  any.  Thereafter,  go  to  step  2. 

Step  2 :  Select  a  suitable  deployment  for  the  sensor  mechanisms 
contained  in  the  combination  obtained  from  step  1  above. 
Thereafter,  go  to  step  3. 

Step  3:  Compute  the  (objective  function)  value  for  the  deployment 
of  step  2 .  If  that  value  is  infinity  (which  means  that  the 
deployment  does  not  cover  the  space) ,  then  go  back  to  step  2  for 
another  deployment,  unless  "chances”  are  exhausted.  (The  number  of 
"chances"  is  set  at  heuristic  initialization) .  Thereafter,  go  to 
step  4. 

Step  4 :  Do  updating  and  book-keeping  to  take  into  account  the  value 
obtained  in  step  3.  Check  stopping  condition .  If  stopping 
condition  is  satisfied,  then  STOP;  otherwise,  go  back  to  step  1, 
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with  information  to  guide  the  next  "skip  and  search". 

The  description  above  gets  formalized  as  the  underlined  basic 
operations  are  explained  in  details.  We  will  now  describe  each  of 
these  operations. 

(a)  "Systematically  Skip  and  Search" 

This  is  essentially  a  line  search,  except  that  the  set  over 
which  this  is  done  is  the  discrete  set  {1,  .  .  .,2“-!}  coded  in  binary 
numbers.  This  set  is  constructed  as  follows.  First,  all 
combinations  of  the  given  sensor  mechanisms  (m  of  them)  are  "lined 
up"  (denumerated)  using  a  binary  coding.  For  example,  let  us 
suppose  m=3 .  Then  the  binary  code  line-up  of  all  combinations  of 
the  three  sensor  mechanisms  is  given  by  the  following: 

001  010  011  100  101  110  111 

I  i  I  I  I  I  I 

{1}  {2}  {1,2}  {3}  {1,3}  {2,3}  {1,2,3} 

Note  that  the  binary  numbers  above  form  an  increasing  sequence  for 
which  the  n-th  term  is  obtained  by  adding  1  to  the  {n-l)-th  term; 
in  decimal  numbers,  that  sequence  is  1,  2,  ...  7. 

The  rule  for  assigning  binary  numbers  to  combinations  of 
sensor  mechanisms  is  as  follows: 

Rule  3.1:  Suppose  a  combination  of  sensor  mechanisms  is  given.  For 
i=l...m,  if  the  i-th  sensor  mechanism  is  included  in  the  given 
combination,  then,  counting  from  right  to  left,  set  the  i-th  digit 
(of  the  binary  nximber)  to  1;  otherwise  set  it  to  0.  {The  inverse  of 
this  operation  is  obvious) . 

Using  Rule  3.1  above,  all  combinations  of  sensor  mechanisms  are 
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lined  up  in  binary  codes  between  1  and  2“-l.  Bence  a  search 
through  all  combinations  of  sensor  mechanisms  becomes  a  search  over 
the  set  {1,  .  .  .  ,2’-l} . 

But  this  "skip  and  search"  does  not  cover  all  of  2“-l} 
every  time.  The  search  confines  itself  to  k-subsets  (of  the  set  of 
sensor  mechanisms)  only.  After  completing  the  search  over  k- 
stibsets,  it  goes  on  to  (k+1) -subsets,  and  thereafter  goes  on  to 
(k+2) -subsets,  and  so  on.  At  the  initialization  of  the  heuristic, 
the  user  is  asked  to  provide  a  starting  k,  if  he  can;  otherwise  m/2 
or  (m+l)/2  is  used  as  initial  k.  This  initial  k  is  the  means  by 
which  the  heuristic  accepts  initial  solution  guesses  from  the  user . 

Details  of  how  this  discrete  line  search  starts  off,  and  how  it 
skips  around  (inside  the  set  of  binary  coded  integers  described 
above)  to  find  the  next  siibset  of  sensor  mechanisms  will  be  given 
later  in  this  section.  For  now,  we  assume  that  it  does  find  the 
next  subset,  that  is,  a  combination  of  sensor  mechanisms. 

(b)  "Select  a  Suitable  Deployment" 

This  is  a  random  selection.  Suppose  d^  is  the  total  nximber  of 
deployment  options  for  the  i-th  sensor  mechanism.  A  random  positive 
integral  nmnber,  say  r,  is  generated,  and  then  deployment  option 
(r  mod  di)+l  is  chosen  for  the  i-th  sensor  mechanism. 

This  way  of  selecting  deployments  has  performed  fairly  well  in 
our  computational  experiments .  This  aspect  of  our  heuristic  is 
amenable  to  certain  advanced  results  in  probability  and  measure 
theory,  but  that  is  not  of  much  relevance  to  our  immediate  purpose 
here. 
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(c)  "Compu'te  the  Objective  Function  Value" 

This  reduces  to  solving  a  standard  set-covering  problem.  Any 
reasonable  heuristic  for  the  standard  set-covering  problem  could  be 
used  here.  We  chose  a  greedy  heuristic,  but  we  know  that  a  divide- 
and-conquer  method  should  be  preferable  on  parallel  computers . 
Again,  computational  experiments  indicate  that  this  is  a  good 
choice . 

Now,  the  given  deployment  of  sensor  mechanisms  may  not  give  a 
cover  for  the  space.  If  the  given  deployment  covers  the  space,  then 
a  finite  value  is  returned  by  the  greedy  heuristic;  otherwise  the 
value  of  this  combination  is  set  to  infinity.  In  the  case  that  the 
given  deployment  does  not  give  a  cover,  the  corresponding 
combination  of  sensor  mechanisms  may  get  another  chance.  The 
deployment  generation  process  described  above  may  be  repeated  for 
the  combination,  unless  the  ntamber  of  "chances"  for  the  combination 
has  reduced  to  0.  The  number  of  chances  is  determined  at  heuristic 
initialization;  it  is  the  maximum  ntunber  of  random  deployments  that 
may  be  generated  for  each  combination  of  sensor  mechanisms  each 
time  it  (the  combination)  is  presented  by  the  "skip  and  search" 
step.  For  every  deployment  generated,  the  number  of  chances  reduces 
by  1. 

(d)  "Do  Updating  and  Book-keeping" 

First,  some  terminology.  As  the  integers  1,  2,  ...  ,  2“-l,  are 
points  of  the  real  line,  we  shall  refer  to  their  binary  equivalents 
0...01,  0...10,  ...  ,  1...11,  as  "points".  Recall  that  each 
of  these  points  corresponds  to  a  subset  of  the  given  set  of  sensor 
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mechanisms .  For  a  point  p,  the  objective  function  value  shall  be 
denoted  by  f (p) . 

Our  heuristic  always  maintains  a  3-point  v-pattern  —  points 

Pl^  Ph/  Pr  such  that  Pi,  <  Pm  <  Pr  aud  f  (Pl)  >  ^  (Pk)  <  ^  (Pr)  • 

Figure  3.1  illustrates  the  3-point  v-pattern.  The  combination  of 
sensor  mechanisms  from  "skip  and  search"  step  corresponds  to  some 
point,  say  p^,  inside  the  interval  [Pl/Pr]  •  "Updating"  consists  of 


Figure  3.1  —  3 -point  v-pattern 

replacing  either  pj,  or  Pr  with  p^  in  the  3-point  v-pattem.  Which 
one  gets  displaced  depends  on  the  magnitude  of  ^  relative  to  that 
of  Pi,  and  the  magnitude  of  f  (Pr)  relative  to  that  of  f  (r,)  .  If 
either  (Pt,  <  Pm  «  f  (Rj)  <  f  (Pm)  )  or  (Pi,  >  r,  &  f  (Ph)  >  f  (R() ) .  then 
Pr  is  displaced;  otherwise  p^  is  displaced. 

"Book-keeping"  consists  of  recording  the  points  constituting 
the  current  3-point  v-pattern  and  their  objective  function  values. 
It  also  includes  recording  information  about  the  "best"  point  so 
far;  this  is  necessitated  by  the  heuristic's  need  to  start  off  and 
"restart"  in  a  certain  manner  (see  (f)  below)  .  As  it  goes  from 
k-sxibsets  to  (k+1) -subsets,  the  heuristic  uses  a  restart  similar  to 
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the  first  heuristic  initialization  (see  (f)  below) . 

(e)  "Check  Stopping  Condition" 

Assume  the  last  "skip  and  search"  step  consists  of  searching 
among  k-subsets .  If  k=m,  our  heuristic  computes  the  objective 
function  value  for  the  point  (1,.../1),  prints  all  results,  and 
STOPS.  Suppose  k  <  m.  If  current  p^^  is  the  only  k-subset  point  in 
the  interval  [pj,,  Pr]  /  then  a  switch  to  (k+1) -stabsets  is  indicated 
as  the  heuristic  returns  to  step  1;  otherwise  the  heuristic  returns 
to  step  1  with  information  that  the  search  is  still  on  k-subsets . 

(f)  More  details  on  the  Line  Search 

Let  L(k)  denote  the  binary  number  with  I's  as  the  first  k 
digits  and  O's  as  the  last  m-k  digits,  counting  from  right  to  left, 
and  R(k)  the  binary  number  with  O's  as  the  first  m-k  digits  and  I's 
as  the  last  k  digits.  Lemma  3.1  below  indicates  how  we  will  use 
L(k)  and  R(k)  . 

Definition  For  the  binary  code  sequence 
0...01,  0.,.10,  0...11,  .  .  .  ,  1...11, 
a  term  shall  be  called  a  k-term  if  it  has  exactly  k  I's  in  its 
digits . 

Lemma  3.1  All  k-terms  are  contained  in  the  interval  [L(k),R(k)]. 
Proof:  This  becomes  obvious  when  one  remembers  that  the  n-th  term 
is  obtained  from  the  (n-l)-th  term  by  adding  1.  ## 

To  start  searching  through  k-subsets,  our  heuristic  sets 
Pi.=L(k),  PR=R(k),  and  pj,  equal  to  a  k-term  close  to  (L(k) +R{k) ) /2 . 
Lemma  3.2  below  shows  how  pj,  may  be  obtained. 

Lemma  3 . 2  (Skipping  Lemma)  Suppose  k  <  m. 
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(i)  Let  p(k+l)  be  a  (k+l)-term.  Obtain  a  k-term  q(k)  from  p(k+l) 
by  setting  to  0  the  rightmost  nonzero  digit  of  p(k+l) .  Then,  there 
is  no  other  k-term  between  q(k)  and  p(k+l) . 

(ii)  Let  u(k)  be  a  k-term.  Obtain  a  (k+l)-term  w(k+l)  from  u(k)  by 
setting  to  1  the  rightmost  zero  digit  of  u(k).  Then,  there  is  no 
other  (k+l)-term  between  u(k)  and  w(k+l) . 

Proof:  Express  terms  in  decimal  niimbers,  and  the  desired 
conclusions  follow  from  straightforward  comparisons .  ## 

When  it  switches  from  k-subsets  to  (k+1) -subsets,  the  heuristic 
restarts  the  search  the  same  way  it  starts  the  search  through  k- 
subsets,  except  that  k  is  replaced  by  k+1. 

When  the  heuristic  starts  or  restarts,  the  objective  function 
values  may  not  give  the  desired  v-pattern.  This  apparent  difficulty 
is  overcome  by  saving  the  best  of  the  objective  function  values,  if 
it ' is  good  enough  (see  "Book-keeping"  under  (d) ) ,  and  thereafter 
setting  f  (Pl)  and  f  (Pr)  artificially  so  as  to  obtain  a  v-pattern. 

After  its  initialization,  the  heuristic  applies  lemma  3.2  to 
the  integral  part  of  (Pl+Pm) /2  or  (Pm+Pr)/2  to  obtain  p^.  It  is  this 
application  of  lemma  3.2  that  constitutes  "skipping"  in  this 
heuristic.  The  kind  of  skipping  done  here  is  similar  to  what  might 
be  done  by  a  procedure  doing  discrete  line  search  on  the  set  of 
integers  divisible  by  7,  then  on  the  set  of  integers  divisible  by 
9,  thereafter  on  the  set  of  integers  divisible  by  11,  and  so  on. 

To  conclude  this  description  of  our  heuristic,  we  now  s'cate  its 
computational  complexity  (that  is,  its  speed) .  First,  our  earlier 
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desciripl^ion  of  "Updating"  needs  a  little  adjustment  so  as  to  ensure 
a  certain  speed-up.  Recall  that  p^  is  obtained  from  (Pi.+Pm)/2  or 
(Pm+Pr)  /2 1  and  then  a  new  v-pattern  is  formed .  Let  us  refer  to  the 
integral  part  of  (Pl+Pm) /2  or  (Pk+Pr)/2  as  the  "source"  of  p^.  To 
ensure  that  it  discards  about  1/4  of  current  interval  with  every 
"skip  and  search"  step,  the  heuristic  must  replace  p^  with  its 
"source"  as  it  (that  is,  p^)  is  about  to  become  a  part  of  the  3- 
point  pattern,  and  regard  f  (p^)  as  the  function  value  for  this 
"source".  This  adjustment  ensures  that  the  heuristic  runs  in 
o  ( I S  |m^logi,32,)  time,  since  our  greedy  procedure  for  set-covering 
runs  in  o(|S|m)  time.  Here,  |S1  is  the  total  niunber  of  vertices  of 
the  cubical  subdivision  of  S  given  in  our  first  technical  report. 

4 .  Corresponding  Computer  Program 

The  major  subroutines  of  our  computer  program  are  close  to  the 
heuristic  operations  described  in  section  2 .  The  flow  charts  of 
Figures  4 . 1  ana  4 . 2  display  the  logic  of  this  computer  program,  and 
the  following  table  describes  a  useful  correspondence  between  the 
heuristic  operations  and  the  program  subroutines. 

Two  simplifying  assumptions  have  been  made.  First,  we  assume 
that  the  space  S  is  a  rectangular  volume  in  R^.  We  also  assume  that 
each  surveillance  sensor  mechanism  has  a  spherical  coverage  or 
span.  These  two  assumptions  are  essentially  assumptions  about 
sensor  technologies,  and  could  be  discarded  without  adversely 
affecting  the  basic  features  of  the  computer  program.  Such 
technological  assoimptions  will  be  thoroughly  considered  in  our 
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follow-up  research  work. 

Heurishic  Operation  Program  Subroutine 

setup  0 

skip  and  search  < - >  restart  () 

nextpt ( ) 

compute  objective  function  < - >  ptvalue() 

updating  and  book-keeping  < - >  update {) 

Computational  experiments  with  this  program  have  been 
encouraging.  The  results  of  these  computational  eaqperiments  shall 
be  fully  described  in  our  next  technical  report. 
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Q,  4,1  -  Flow  Chart  for  MainC  ) 
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/*  discrete  line  search  heuristic  main  driver*/ 
fdefine  CHANCES  1 
tdefine  BIG  5000 
fdefine  LIM  20 
fdefine  RIM  80 
finclude  "malloc.h" 
finclude  "stdlib.h" 

'  finclude  "stdio.h" 
void  setup(void); 
void  nextpt ( void ) ; 
void  ptvalue ( void ) ; 
void  update ( void ) ; 
void  restart { void )  ; 
void  output( void) ; 
void  binen ( int ) ; 
int  probe(int,  int.  int,  int); 
int  space [ 5 ] ; 

float  sensor[ 2 ] [LIM] ,  deploy [ 5 ] [RIM] ; 

int  subset[LIM],  index[LIM]; 

int  ns,  kstrt,  jj,  objflag; 

int  pnt[3],  newpt,  gap; 

float  value[3],  newval,  tempv; 

float  resultl [ 5 ] [LIM ] ,  aresult [ 5 ] [LIM ] ,  cresul t [ 5 ] [ LIM ] ; 

int  1,  w,  h,  pieces,  seed; 

char  fname[6],  fnamel[7],  fname2[7]; 

main  (  ) 

{ 

FILE  *fptr,  *fptrl,  *fptr2; 
int  i; 

printf  (  " _  _ _ ^\n\n"  )  ; 

print£{"The  set  of  data  for  this  procedure  exists  as  a  DATA\n"); 
printf { "FILE .  You  are  required  to  give  the  NAME  of  this  file.\n\n"); 
printf ("You  may  not  use  more  than  6  characters  for  this  file\n"){ 
printf ( "name . \n" ) ; 

printf  (  " _ _ _ _ _ ^\n\n\n"  )  ; 

printf ( "Please  give  the  FILE  NAME  -  no  more  than  6  characters  -->  "); 
f flush ( stdin ) ;  scanf("%s",  fname);  printf (" \n\n" ) ; 
fnamel [0 ] » ' s ' ;  fname2 [0 ]  =  ' d ' ; 
for  (j-0;  j<6;  j++) 

{ 

fname 1 [ j+1 ] -fname [ j ] ;  fname 2 [ j+1 ] -fname [ j ] ; 

} 

fptr-f open ( fname ,  "rb"); 

fread{spaca,  sizaof ( space ) ,  1,  fptr);  f close { fptr ) ; 
f ptrl=f open ( f namel ,  "rb"); 

f read ( sensor ,  sizeof ( sensor ) ,  1,  fptrl);  f close ( fptrl ) ; 
f ptr2=f open ( f name2 ,  "rb"); 

f read ( deploy ,  sizeof { deploy ) ,  1,  fptr2);  f close { fptr2 ) ; 
setup! ) ; 

for  (  jj  »  kstrt;  jj  <  ns;  jj++  ) 

{  nextpt ( ) ; 

while  {  gap  !-  0  ) 

{  ptvalue();  update();  nextpt();  } 
restart( ) ; 

} 

output ( ) ; 

} 


- -  f 

/*  initiates  the  line  search  iterations  */ 
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void  setup(vold) 

{ 

int  i.  j.  k.  temp,  tempi,  sum; 
l»space [0 ] /space [ 3 ]  +  1; 
w-space[l]/space[3i  +  1; 

h-space [ 2 ] /space [ 3 ]  +  1;  pieces-l*w*h;  seed-1; 
ns-space[4l :  sum-0: 

for  (j-0;  j<ns;  j++) 

{  index[ j ]»sum;  temp-sensor [ 1 ][ j 1 ;  sum+»temp;  } 
printf ( "Please  SELECT  an  objective  for  this  surveillance  problem\n\n" ) ; 
printf("  1  --  To  Minimise  Maintenance  Costs\n\n”); 

printf ("  2  --  To  Minimise  Maintenance  Task  Time\n\n"): 

printf ("  YOUR  SELECTION  -->  "); 

scanf("%d",  &objflag);  printf (" \n\n" ) ; 

printf ( "Would  you  guess  the  number  of  surveillance  sensor  mechanismsNn" ) 

printf ("that  may  suffice  in  covering  the  given  space?  If  so,  ENTERNn"); 

printf ("the  number  here,  OTHERWISE  ENTER  the  number  zero  (0)  -->  "); 

scanf("%d",  &k);  printf (" \n\n" ) ; 

if  (k--0  II  k>ns-2)  kstrt-ns/2;  else  kstrt-k; 

jj=kstrt;  temp-0; 

for  (  j-0;  j<kstrt;  j++)  /‘first  left-point  and  corresponding  subset*/ 

{  subset[jl-l;  terap=2* temp+1 ;  } 
for  (  j=kstrt;  j<ns;  j++)  subset[jl-0; 

ptvalue();  value [0 ]=newval ;  pnt(0]=temp;  tempv-newval ; 
for  (i-0;  i<5;  i++) 

for  (j=0;  j<ns;  j++)  aresult [ i 1 [ j ] -cresul t [ i ] [ j ] ; 
for  (  j-0;  j<(ns-kstrt) ;  j++)  subset[j]=0;  /*  first  right-point  and*/ 
for  (  j-(ns-kstrt) ;  j<ns;  j++)  subset[j]=l;  /‘corresponding  subset*/ 
for  {  j=kstrt;  j<ns;  j++)  temp  *«  2; 
ptvalue();  valueI2]-newval;  pnt[2J=temp; 
if  (newvai<tempv) 

{  tempv-newval; 

for  (i-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  aresult [ i J [ j ] -cresult [ i ][ j ] ; 

} 

pnt[l]=(pnt[01+pnt[2l ) ;  binen(pnt(ll ) ; 
ptvalueH;  value [  1  j-newval ; 
for  (i-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  resultl [ i ] [ j ]-cresult ( i  ]  [ j  ] ; 

} 


- 

/*  computes  the  next  point  and  corresponding  subset  */ 
void  nextpt(void) 

{  int  p0,  pi,  p2; 
p0»pnt[0];  pl=pnt[l];  p2=pnt[2]; 
if  (pl-p0  >-  p2-pl) 

{  if  (pi  >  p0  +  1)  {  newpt  »  (pl+p0)/2;  gap-1;  } 
else  gap  »  0; 

} 


if  (pl-p0  <  p2-pl) 

{  if  (p2  >  pi  +  1)  {  newpt  »  (pl+p2)/2;  gap-2;  } 
else  gap-0; 


} 


binen ( newpt )  ; 


} 


- 

/*  updates  points  and  their  values  */ 
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void  update(vold) 

{  int  i.  j; 

if  (newval  <  value[l]  ) 

{  for  (i-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  resultl [ i] [ j ]-cresult [ i I [ j ] ; 
if  (gap«l) 

{  value [2 J-valUG [ 1 ] ;  value [IJ-newval; 
pnt[2]«pnt[l] ;  pnt [ 1 ]»newpt; 

} 

if  (gap=»-2) 

{  value [0 ] »value [ 1 ] ;  value [ 1 J =newval ; 
pnt[0]»pnt[l] ;  pnt[ 1 ]=newpt; 

} 

} 

else 

{  if  (gap--!) 

{  value[0]*newval;  pnt[0]=newpt;  } 
if  (gap»»2) 

{  value [ 2 ]*newval ;  pnt[2]»newpt;  } 

} 

} 


- 

/*  does  a  binary  enumeration  */ 
void  binen{int  p) 

{  int  XX,  j; 

XX-p; 

for  (j=0;  j<ns;  j++)  subset[j]=0; 
j*0;  - 

while  {  XX  >  1  fieSt  j  <  ns  ) 

{  subset[j]  ■  xx\2i 

XX  «  (XX- (xx%2) ) /2 ;  j+-ls 

} 

subset [ j ] »1 ; 


/ . . . . */ 

/*  restarts  the  line  search  iterations  */ 
void  restart ( void) 

{ 

int  i,  j,  k,  p,  pi,  temp,  flag; 
if  ( value[l ]<tempv) 

{  tempv*value [ 1 ] ; 
for  (i*0;  i<5;  i++) 

for  (j»0;  j<ns;  j->-+)  aresult  (  i  ][  j  ] -resultl  [  i  ][  j  ] ; 

} 

p*0 ; 

if  (jj=>*ns-l) 

{  for  (j-0;  j<ns;  j++) 

{  subset[jl-l;  p»(2*p)  +  1  ;  } 
ptvalue ( )  ; 
if  (newval  <  tempv) 

{  tempv  »  newval; 
for  (i*0;  i<5;  i++) 

for  {j-0:  j<ns;  j++)  aresult [ i ] [ j ] -cresult [ i  ]  [ j  ]  ; 

} 

} 

if  (jj<ns-l) 

{  P-0; 

for  (  j-0;  j<jj+l;  j++)  /*  left-point  and  corresponding 


subset*/ 
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{  subset[j]»l;  P-2*P+1;  } 
for  (  j-jj+1;  J<ns;  j++)  subset[j]-0; 
ptvalue();  value[0]«newval;  pnt[01-p;  tempv-newval ; 
for  (i-0;  i<5;  i++) 

for  {j"0;  j<nsj  j++)  aresult [ i] [ j ]-cresult [ i 1 [ j ] ; 
for  (  j-0;  j< ( ns- ( j j+1 ) ) ;  j++)  subsetlj]»0;  /*  right-point  and  */ 
for  (  j-(ns-( j j+l) ) ;  j<ns;  j++)  subset[j]»l;  /‘corresponding  subset*/ 

(  j-jj+1;  j<ns;  j++)  p  *-  2; 
ptvalue():  value t 2 ]-newval ;  pnt(2]-p; 
if  (newval<teinpv) 

{  tempv-newval ; 

for  (i-0;  i<5;  i++) 

for  {j-0;  j<ns;  j++)  aresult(i]  [  j  ]=*cresult[i]  [  j  ]  ; 

} 

pnt[l ]-{pntE0]+pntl2] ) /2;  binen (pnt [ 1 ] ) ; 
ptvalue{);  value [ 1 ]»newval ; 
for  (i-0;  i<5;  i++i 

for  (j-0;  j<ns;  j++)  resultl [ i 1 [ j ]»cresult [ i ] [ j ] ; 

} 

} 

- 

/*  prints  the  output  on  the  screen  */ 
void  output(void) 

{  int  i,  j; 
char  fn[4]; 
if  (objflag--!) 

{  fn[01-'c';  fn[ll-'0';  fn[2]-'S';  £n[31-'T';  } 
if  (objflag— 2) 

{  fnl0]-'T';  fnll]-'l';  £nI2]-'M';  fn[3]-'E';  } 
if  (  tempv  <  value[l]  ) 

{  value [ 1 ] -tempv; 
for  (i-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++}  resultl [ i ] [ j J-aresult I i ] I j ] ; 

} 

printf("  OUTPUT  STATEMENTNn" ) ; 

printf("  _ \n\n\n"); 

printf("  Input  FILE'S  name  is  ’\s'  \n\n\n’',  fname); 

printf("  Corresponding  to  ic%c%c%c  value  of  %5.2f,  the  f ollowingNn" , 
fn[0],£n[l],  fn[2],  fn[3],  valuell]); 
printf("  table  gives  details  of  recommended  surveillance  sensorSn"); 
printf("  mechanism  deployments\n\n" ) ; 

printf{"  _ _  \n\n"); 

printf ( "Sensor\tRadius\t\tLocation\t\tCost\t\tTime\n\n" ) ; 
for  (j-0;  j<ns;  j++) 

{  if  (  resultl [4 ][ j 1  >  -1  ) 

printfC  %d\t%5.2f\t  (15. 2f.  15. 2f,  15.2f)\t  %4.2f\t  %5.2f\n". 

j  +  1,  sensor [0  ][ j  ] ,  resultl [0 ][ j ] ,  resultl [ 1  ][ j  ] ,  resultl [ 2 ][ j  ]  , 
resultl[3]( j],  resultl(4][j]) ; 

} 

} 

/* . . - . - . / 

/•  checks  whether  the  given  space  may  be  covered  by  a  given  */ 

/*  set  of  sensor  mechanisms,  and,  if  so,  computes  the  'value'  */ 

/*  of  the  cover  relative  to  given  objective  function  */ 

void  select(void) ; 

void  sort ( void ) ; 

void  xyz ( void ) ; 

void  check{void); 
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int  p.  quit; 

Int  polnter[LIM] .  slct[LIM].  part[LIM]; 
float  X.  y.  z.  costime[LIM] : 
void  ptvalue{vold} 

{  int  j,  count;  t*  selects  deployment  */ 

count«0; 

while  (count  <  CHANCES) 

{  for  (j-0;  j<ns;  j++) 

{  cresult[01 [ j 1-0;  cresult ( 1 J ( j ]»0;  cresult [ 2 ] [ j ]»0 ; 
cresultisj [ j i»0;  cresult ( 4 j i*-2 ;  part[jl»0; 

} 

select();  sort{);  p-1;  quit»0;  newval«0; 
while  {  p  <  pieces+1  &&  quit  *-  0  ) 

{  xyz{ )  ;  check( ) ; 

if  (quit  0)  p+*l; 

} 

if  (quit  ==  0)  count-CHANCES+1 ;  /*  a  cover  found  */ 

else  count+*l; 

} 

if  (quit-»l)  /*  not  a  cover  */ 

{  newval=BIG; 

for  (j=0;  j<ns;  j++) 

{  cresult[0] [ j ]=0;  cresultl 1 ] [ j J=0;  cresult [ 2 ][ j ] =0 ; 
cresult[3i[ji=0;  cresult[4j[ji=>-2; 

} 

} 

} 


/* - - - - - - - - - V 

/*  selects  a  deployment  of  a  given  set  of  sensor  mechanisms  */ 
void  select(void) 

{  int  j,  r,  num; 

for  (j=0;  j<ns;  j++) 

{  num»seed*7;  r-sensor [ 1 ] [ j ] ; 

slct ( j ]» ( subset [ j ] ) * ( num%r ) ;  seed+»l; 

} 

} 

- 

/*  sorts  deployed  sensor  mechanisms  by  costs  or  by  time  •/ 
void  sort(void) 

{  int  i,  j,  k,  col,  tempi,  count; 
float  temp; 

if  (objflag=»l)  k-3;  else  k»4; 

/*  copying  into  costime  */  i»0; 
for  (j»0;  j<ns;  j++) 

{  if  { subset [ j ] ==1 ) 

{  col»index[ j ]+slct[ j ] ; 

costime [ i j “deploy [ k] [col J ; 
pointer[ii*j :  i+=l; 

} 

} 

/•  sorting  */  count*jj; 
while  (count>0) 

{  for  (i*0;  i<count;  i++) 

{  if  ( costime [ i ] >costime [ i  +  1 J ) 

{  temp“costime ( i ] ;  costime [ i ]*costime [ i+1 ] ; 
costime [i+1]* temp;  templ=pointer [ i ] ; 
pointer [ i ] “pointer] i+1 1 ;  pointer] i+1 ] *  tempi; 

} 
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> 

count--l ; 

> 


. . */ 

/*  computes  the  origin  of  the  piece  corresponding  to  p  */ 
void  xyz(void) 

{ 

int  num,  a,  b,  c.  temp; 
num»p-l;  temp  «  l*w; 
c  ■  num/temp;  num  temp; 
b  -  num/1;  a  -  num%l; 

temp=space [ 3 ] ;  x»a*temp;  y=b*temp;  2=c*temp; 

} 


/*  . / 

/*  verifies  that  the  point  (x.y.z)  is  covered  */ 
void  check(void) 

{  float  a,  b,  c,  r,  sa,  sb,  sc,  sr; 
int  i,  j,  col,  flag; 
flaga0;  i»0; 
while  {flag==0  &&  i<jj) 

{  jspointer [ i ] ;  col=index[j]  +  slct[j]; 
r= sensor (0 ][ j ] ;  a=deploy [0 ] [col ] ; 
b*deploy[l j [col] ;  c=deploy [ 2 ] [col ] ; 

sr=r*r;  sa=(a-x)  *  (a-x) ;  sb» (b-y) •  {b-y ) ;  sc=«(c-2)  *  (c-z) ; 
if  ((sa+sb+sc)  <=  sr) 

{  flag=-l;  cresult[0]  [  j]=a;  cresult[l]  [  j  ]-b; 
cresultt2I[ j]=c;  cresult [3] [ j l-deploy[3 J [col ] ; 
cresult[4] [ j ]-deploy[4l[colI; 
if  (part[j]— 0) 

{  if  (objflag=-l)  newval+-cresult[3] [ j ] ; 
if  (objflag»-2  &&  newval<cresult[4 ] [ j ] ) 
newval-cresult [4 ] [ j ] ; 
part[ j ]-l; 

} 

} 

i+-l; 

} 

if  (flag-»0)  quit-1; 
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SOURCE  CODE  FOR  COMPUTER 

PROGRAM 


23 


I 


/*  A  GENERALISED  SET-COVERING  PROCEDURE  FOR  SURVEILLANCE  MAINTENANCE  */ 
void  datafls( void) ; 
void  searchl ( void) ; 
void  main (void) 

{ 

char  choice; 
menui 

prlntf("\n\n  MAIN  MENU  \n\n"); 

prlntfi"  _  \n\n\n\n"); 

prlnt£("  SELECT  \n\n\n"); 

prlntf("  A  --  To  Create  or  View  a  Data  Flle\n\n"); 

prlntfi"  B  --  To  Run  this  Space-covering  Procedure\n" ) ; 

prlntf("  on  an  Existing  Data  Flle\n\n"); 

prlntf("  C  --  To  Quit  and  Return  to  System\n\n\n" ) ; 

prlntf("  YOUR  SELECTION  -->  ");  cholce-getche ( ) ; 

prlnt£("\n\n") ; 
swltch(  choice  ) 

{ 

case  'a' I 
case  ' A' i 

data£ls():  prlntf ( "XnXnPress  ENTER  to  contlnue\n\n\n" ) ; 
getch();  goto  menu; 
case  'b't 
case  'B ' I 

searchl();  prlntf { " \n\nPress  ENTER  to  contlnue\n\n\n" ) ; 
getch();  goto  menu; 
case  'c ' I 
case  'C'l 
exlt(0 ) ; 
defaults 
goto  menu; 

} 

} 

/*  ee§e§§eeeee§§§eeeee§ee§e§eee§ee0@§eeeee§§ee§e§§ee@§e@§@@  */ 

void  w_data ( void ) ; 
void  v_data( void) ; 
void  datafls{ void) 

{ 

char  choice; 
menui t 

prlntf ("\n\n  MENU  FOR  DATA  FILESNn"); 

prlntf  (  " _ \n\n"  ) ; 

prlntf ("  SELECT\n\n" ) ; 

prlntf ("  1  --  To  Create  a  New  Flle\n\n"); 

prlntf ("  2  --  To  View  an  Existing  Flle\n\n"); 

prlntf ("  3  --  To  Return  to  MAIN  Menu\n\n"); 

prlntf ("  YOUR  SELECTION  -->  ");  cholce-getche () ; 

prlntf { "\n\n" ) ; 
switch!  choice  ) 

{ 

case  'I't 

w_data{);  prlntf ( "\n\nPress  ENTER  to  contlnue\n\n\n" ) ; 
getch();  goto  menui; 
case  '2't 

v_data();  prlntf { "\n\nPress  ENTER  to  contlnue\n\n\n" ) ; 
getch();  goto  menui; 
case  '3' I 
return; 
defaults 


z 

goto  nenul; 

} 

} 

/*0§§e§§§§ee§ee§eeee@§eeee§§§e@ee@§e§§eeeeeee§e@e0ee§§e§e§  */ 

tdeflne  LIM  20 
tdefine  RIM  80 
ilnclude  "malloc.h" 
tinclude  "stdllb.h" 
ilnclude  "stdio.h" 

Int  space [5]; 
float  sensor[2 ] [LIM] ; 
float  deploy [ 5 ] [RIM] ; 
void  spc ( void)  ; 
void  snsr(vold); 
void  dply(vold); 

FILE  *fptr,  *fptrl,  *fptr25 

char  fnaine[6].  fnamel[7],  fname2[7]: 

void  w_data(vold) 

{ 

Int  j; 

prlntf  (  " _ _ _ \n\n"  ) ; 

prlntf{"The  set  of  data  that  you  are  about  to  enter  will  formXn"); 
prlntf ("a  FILE.  You  are  required  to  give  a  NAME  to  this  flle.\n\n"); 
prlntf ("You  may  not  use  more  than  6  characters  for  the  f lle\n" ) ; 
prlntf ( "name . \n" ) ; 

prlntf  (  " _ \n\n\n"  ) ; 

prlntf { "Please  give  the  FILE  NAME  -  no  more  than  6  characters  -->  "); 
f f lush( stdln) }  scanf("%s",  fname);  prlntf (" \n\n" ) ; 
fnamel [0]- ' s ' j  fname2 [0 ]- ' d ' ; 
for  (j-0;  j<6;  j++) 

{ 

fnamel [ j  +  1 ] -fname [ j ] ;  fname 2 [ j-M ] -fname [ j ] 5 

] 

spc();  snsr();  dply(); 

} 

/*  Information  on  space  under  surveillance*/ 
void  spc (void) 

{ 

Int  j,  k; 

prlntf  (  "  _ \n\n"  ) ; 

prlntf ("  You  will  be  asked  to  provide  Information  about  the\n"); 
prlntf ("  space.  S,  under  surveillance,  and  the  potential  Intruder ' s\n" ) 
prlntf ("  size.  We  shall  assume  that  S  Is  a  rectangular  volume  \n"); 
prlntf("  located  at  the  origin. (0.0.0) .  of  the  3-dlmenslonal  space. \n") 

prlntf  (  " _ \n\n"  ) ; 

prlntf ( "Please  give  LENGTH.  WIDTH.  HEIGHT  of  S.  Irw.h  -->  "); 
fflush( stdln ) ; 

scanf ( "%di %d t %d” .  &space[0],  &space[l].  &space[2])i 
prlntf ( "\n\n" ) ; 

prlntf ( "Please  give  SIZE  OF  POTENTIAL  INTRUDER  -->  "); 
fflush(stdln) ; 

scanf("ld",  &space[3]):  prlntf ( "\n\n" ) ; 

prlntf ("Please  give  TOTAL  NUMBER  OF  AVAILABLE  \n"); 

prlntf ("SURVEILLANCE  SENSOR  MECHANISMS  -->  ");  f flush ( Stdln ) ; 

scanf("%d",  &space[4])i  prlntf {" \n\n\n\n" ) ; 

If  (  (fptr  -  fopenj fname.  "wb" ) ) --NULL) 

{  prlntf ( "Can ' t  open  %s  n".  fname);  exlt();  } 
fHrlte(8pace.  slzeof (space  .  1.  fptr); 
fclose ( fptr) ; 


} 

/*  DATA  ON  SENSORS*  SPAN  RADIUS.  AND  DEPLOYMENT  OPTIONS*/ 
void  snsr(void) 

{ 

int  j; 

prlntf  (" _ \n\n")  ; 

printf{"you  will  next  supply  information  about  the  coverage  \n"); 
printf("of  each  sensor  mechanism  —  the  span,  and  the  number  of\n")j 
printf ( "possible  positioning  or  deployment  of  each  sensor  \n"); 
printf ( "mechanism.  Each  sensor  mechanism  'covers'  a  sphereNn"); 
printf ( "whose  centroid  will  play  the  role  of  deployment  location . \n" ) 

printf  (" _ \n\n\n"  ) ; 

for  (j-0;  j<space[4];  j++) 

{ 

printf ( "Please  give  the  COVERAGE  RADIUS  of  sensor  mechanism\n" ) ; 
printf ( "number  %d  -->  ",  j+1);  f flush ( stdin ) ; 
scanf("%f",  fitsensor [0 ] [ j ] ) ;  printf ( "\n\n" ) ; 

printf ( "Please  give  the  NUMBER  OF  POSSIBLE  DEPLOYMENTS  \n" ) ; 
printf ("of  sensor  mechanism  number  %d  — >  ",  j+1);  f flush ( stdin ) ; 
scanf("%f",  fiisensor [ 1 ] [ j ] ) ;  printf ( "\n\n" ) ; 

} 

if  (  (fptrl  -  f  open  ( fnamel ,  "wb"))«NULL) 

{  printf ( "Can ' t  open  %s\n",  fnamel):  exit();  } 
fwrite ( sensor ,  sizeof ( sensor ) ,  1,  fptrl); 
f close ( fptrl ) ; 

} 

/*DATA  ABOUT  DEPLOYMENT  OPTIONS*/ 
void  dply(void) 

{ 

int  j,  Jc,  m; 

printf  (  " _ \n\n"  )  ; 

printf ("You  will  next  supply  information  concerning  deployment\n" ) ; 
printf ( "locations,  deployment  costs,  and  the  amount  of  time  \n"): 
printf { "needed  to  complete  each  deployment  as  a  task,  that\n"); 
printf ("is,  the  amount  of  time  needed  to  move  each  sensorNn"); 
printf ( "mechanism  from  current  location  to  deployment  location. \n" ) ; 

printf  (  " _  \n\n\n"  )  ; 

■»0; 

for  (j-0;  j<space[4];  j++) 

{ 

for  (k-0;  k<sensor[ 1 ] [ j ] ;  k++) 

{ 

printf ( "Please  give  the  DEPLOYMENT  LOCATION,  x,y,z  of  the  \n" ) 
printf ("%d  th  deployment  of  sensor  mechanism  %d  -->  ",  k+1,  j+1); 
fflush(stdln) ; 

scanf ( "%f ,%f ,%f " ,  &deploy[0 J [m] ,  &deploy[ 1 ] [m] ,  Sdeploy [ 2 ] [m J ) ; 
printf ( "\n\n" ) ; 

printf ("Please  give  the  COST  and  DEPLOYMENT  TASK  TIME.  cst.  \n"); 
printf ("of  the  %d  th  deployment  of  sensor  mechanism  %d  -->  ", 
k+1,  j+1);  f flush ( stdin ) ; 

scanf  (  "%f  *  %f"  ,  &deplc  ,'[3  ]  [m] ,  &deploy[4  ]  [m] )  ;  printf  (" \n\n"  ) ; 
m+-l ; 

} 

> 

if  (  (fptr2  -  fopen(  fnam-^  2.  wb"))— NULL) 

{  printf ( "Can' t  open  lie  please");  exit();  } 
fwrlte(deploy.  sizeof (de  Lov  .  1,  fptr2); 
f close ( fptr2 ) ; 

} 

/*  eeee§ieeeee@§eeee§§@ee§@@  i§@eee§eeee§§0§eeeeee§eeeee§e§e§e@e  */ 


void  v_data{void) 

{ 

Int  j,  k.  m; 

prlntf  {  " _ \n\n"  ) ; 

prlnt£("The  set  of  data  that  you  are  about  to  view  exists  asVn*); 
prlntf ("a  FILE.  You  are  required  to  give  the  NAME  of  this  file.\n\n"); 
prlntf ("You  may  not  use  more  than  6  characters  for  this  file\n"); 
prlntf { "name . \n" ) ; 

prlntf  ( " _ \n\n\n"  ) ; 

prlntf ( "Please  give  the  FILE  NAME  -  no  more  than  6  characters  -->  "); 
f flush! stdln) ;  scanf("%s",  fname);  prlntf (" \n\n* ) ; 
fnamel [0 ]- ' s ' :  fname2[0]-'d' ; 

for  {j-0;  j<6;  j++) 

{ 

fnamel [ j+1 ] -fname [ j ] ;  £name2 [ j+1 J-fname [ j 1 ; 

> 

fptr«fopen( fname.  "rb"); 
fread(space.  slzeof ( space ) .  1.  fptr); 
fclose{fptr) ; 

prlntf ("%s  DATA  ABOUT  SPACE,  S,  UNDER  SURVBILLANCBSn" ,  fname); 

prlntf  (  " _ _ \n\n"  ) ; 

prlntf ( "Dimensions!  %d,  %d.  %d\n\n".  space[0].  space[l],  space[2]); 
prlntf !"Size  of  potential  intruder!  %d\n\n\n",  space[3]); 
prlntf ( "Press  any  key  to  continue  file  viewing\n\n\n" ) ; 
getch( ) ; 

prlntf ("%s  DATA  ABOUT  AVAILABLE  SENSOR  MECHANISMS \n " ,  fname); 

prlntf  (  " _ \n\n"  ) ; 

fptrl«fopen( fnamel .  "rb"); 

fread( sensor,  slzeof ( sensor ) ,  1,  fptrl); 

f close ( fptrl ) ; 

prlntf ! "Number  of  sensor  mechanisms!  %d\n\n",  space[4]); 

prlntf ( "sensor  mechanism  radius  of  span  deployment  optlonsVn"); 

for  (j-0;  j<space[4I;  j++) 

prlntf ( "\t%d\t\t%f \t\t%f\n" ,  j  +  1,  sensor [0] [ j ] ,  sensor [ 1  ][ j  ]) : 
prlntf ( "\n\n\n" ) ; 

prlntf! "Press  any  key  to  continue  file  viewing\n\n\n" ) ;  getch(); 
prlntf {"%s  DATA  ABOUT  DEPLOYMENT  OPTIONSNn",  fname); 

prlntf  (" _  Sn\n'')t 

fptr2-fopen( fname2 ,  "rb"); 
fread(deploy,  slzeof (deploy) ,  1,  fptr2); 
f close ( fptr2 ) ;  m-0; 

for  (j-0;  j<spdCe[4]:  j++) 
for  (k-0;  k<sensor[ 1 ] [ j ] ;  k++) 

{ 

prlntf ( "<--DEPLOYMENT  OPTION  %d  of  sensor  mechanism  %d — >\n\n" ,  k+1,  j+1) 
prlntf ( "Location!  (  %f,  it,  %f  )\n\n",  deploy[0] [m] .  deploy! 1 ] 1 m] , 
deploy[2] [m] ) ; 

prlntf ( "Cost!  %f  \t\t",  deploy[3] [m] ) ; 

prlntf ! "Deployment  task  tlmei  %f\n\n\n",  deploy[4] [m] ) ;  m+-l; 

prlntf ( "Press  any  key  to  continue  file  vlewlng\n\n\n" ) ;  getch(); 

> 

} 

/"  e§eeee§0eeee§§§eeee§eeeeeee@0e0eee§e§§eee§eee§e0eeeeeee§§§ee@eee@  */ 

/*  discrete  line  search  heuristic  main  driver*/ 

«deflne  CHANCES  1 
«de£ine  BIG  5000 
«deflne  LIM  20 
ideflne  RIM  80 
ilnclude  "malloc.h" 
ilnclude  "stdllb.h" 


iinclude  "stdlo.h" 
void  setup ( void ) ; 
void  nextpt ( void ) ; 
void  ptvalue(void) i 

void  update(void) :  t*  functions  */ 

void  restart ( void ) ; 
void  output( void) : 
void  binen( int) : 

int  probe(int.  int,  int,  int);  /*  global  variables  •/ 

int  space[5]; 

float  sensor [2 ] [LIM] ,  deploy [5 ] [RIM] ; 

int  subset[LIMi,  index[LIM]; 

int  ns,  kstrt,  jj,  objflag; 

int  pnt[3],  newpt,  gap; 

float  value[3],  newval,  tempv; 

float  resultl [ 5 ] [LIM] ,  aresult [5 ] [LIM] ,  cresult [ 5 ] [LIM] ; 

int  1,  w,  h,  pieces,  seed; 

char  fname[6].  fnamel[7].  fname2[7]; 

void  searchl ( void)  /*  beginning  */ 

{ 

FILE  *fptr,  ‘fptrl,  *fptr2; 
int  j; 

printf  (  " _ \n\n"  )  ; 

printf("The  set  of  data  for  this  procedure  exists  as  a  DATA\n"); 
printf ( "FILE .  You  are  required  to  give  the  NAME  of  this  file.\n\n"); 
printf ("You  may  not  use  more  than  6  characters  for  this  fileVn"); 
printf ( "name . \n" ) ; 

printf  (  " _ \n\n\n"  ) ; 

printf { "Please  give  the  FILE  NAME  -  no  more  than  6  characters  -->  "); 
f f lush( stdln) ;  scanf("%s",  fname);  printf (" \n\n" ) ; 
fnamel [0]« ' s ' ;  fname2[0]>'d' ; 

for  (j-0;  j<6;  j++) 

{ 

fnamel [ j+1 ]-fname [ j ] ;  fname2 [ j+1 ]-fname [ j ] ; 
fptr-f open ( fname ,  "rb"); 

freadCsj'  ce,  slzeof  ( space ) ,  1,  fptr);  fclose ( fptr ) ; 
fptrl«f'  ien(fnamel,  "rb"); 

f read ( stinsor ,  slzeof  ( sensor )  ,  1,  fptrl);  fclose (  fptrl ) ; 
fptr2-fopen(fname2,  "rb"); 

fread(deploy,  slzeof  ( deploy ) ,  1,  fptr2);  f close ( fpti:2 ) ; 
setup! ) : 

for  (  jj  •  kstrt;  jj  <  ns;  jj++  ) 

{  nextpt();  /*  see  flow  chart  */ 

while  !  gap  !-  0  ) 

{  ptvalue();  update();  nextpt();  } 
restart ( ) ; 

} 

output ( ) ; 

> 

- 

/*  initiates  the  line  search  iterations  */ 
void  setup (void) 

{ 

int  i,  j,  k,  temp,  tempi,  sum; 
l-8pace[0]/space[3]  1; 
w-8pace[l]/8pace[3]  *  1; 

h*8pace[2]/space[3]  +  1;  pi  ces-l*w*h;  seed-3*kstrt+l ;  /*  random  seed  * 
n8«space[4] ;  sum-C; 


for  (j-0j  j<ns;  j++) 

{  Index!  j  ]-sum:  temp-sensor [  1  ][  j  ] ;  sum-*— temp;  } 
prlntf ( "Please  SELECT  an  objective  for  this  surveillance  problem\n\n" ) ; 
printf("  1  —  To  Minimise  Maintenance  Costs\n\n"); 

prlntf ("  2  —  To  Minimise  Maintenance  Task  Time\n\n"); 

prlntf {"  YOUR  SELECTION  — >  "); 

scanf("%d",  Stobjflag);  prlntf  (  "\n\n"  ) ; 

prlntf ( "Would  you  guess  the  number  of  surveillance  sensor  mechanlsmsNn" ) 

prlntf ("that  may  suffice  in  covering  the  given  space?  If  so,  ENTERXn" ) ; 

prlntf ("the  number  here,  OTHERWISE  ENTER  the  number  zero  (0)  — >  "); 

scanf("%d",  &k) ;  prlntf ( "\n\n" ) ; 

if  (k--0  II  k>ns-2)  kstrt-ns/2;  else  kstrt-k; 

jj-kstrt;  temp-0: 

for  (  j-0;  j<kstrt;  j++)  /"first  left-point  and  corresponding  subset*/ 

{  subset! j]-l;  temp-2*temp+l;  } 
for  (  j-kstrt;  j<ns;  j++)  subset! j]-0: 
ptvalue():  value!0]-newval:  pnt!01-temp;  tempv-newval; 
for  (1-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  aresult ! 1 ] ! j ]-cresult ! 1 ] ! j ] ; 
for  (  j-0;  j<(ns-kstrt) ;  j++)  subset! j]-0;  /*  first  right-point  and*/ 
for  (  j- (ns-kstrt ) ;  j<ns;  j++)  subset! j]-l:  /"corresponding  subset*/ 
for  (  j-kstrt;  j<ns;  j++)  temp  *-  2; 
ptvalue();  value !2]-newval;  pnt!2]-temp; 
if  (newval<tempv) 

{  tempv-newval; 

for  (1-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  aresult ! 1 ]! j J-cresult! 1 ]! j ] ; 

} 

pnt!  l,]»(pnt!0]+pnt!2] ) /2;  blnen(pnt!l] ) ;  sum-0: 
for  (j*'0;  j<nsi  j++)  sum+-subsetl  j  ] ;  j*0; 

while ( sum< j  j  ) 

{  sum+-l-subset!  j  ] ;  subset!j]-l;  j+-l;  } 
while ( sum> j  j ) 

{  sum--subset ! j ] :  subset! j]-0;  j+»l;  } 
ptvalue();  value ! 1 I-newval; 
for  (1-0;  1<5;  i++) 

for  (j-0;  j<ns;  j++)  resultl !1 ]! j J-cresult !i J! j ] ; 

} 


/*  computes  the  next  point  and  corresponding  subset  */ 
void  nextpt(vold) 

{  int  p0,  pi,  p2,  j,  sum; 
p0-pnt!0];  pl-pntlU;  p2-pntl2]; 
if  (pl-p0  >-  p2-pl) 

{  if  (pi  >  p0  +  1)  {  newpt  -  (pl+p0)/2;  gap-1;  } 
else  gap  -  0; 

} 

if  (pl-p0  <  p2-pl) 

{  If  (p2  >  pi  +  1)  {  newpt  -  (pl+p2)/2;  gap-2;  } 
else  gap-0; 

} 

binen ( newpt ) ;  sum-0 ; 

for  (j-0;  j<ns;  j++)  sum+-subsetl j J ;  j-0; 

while (sum<jj ) 

{  sum-f-l-subset!  j  ] ;  subset!j]-l;  j+-l;  } 
whlle(sum>j j ) 

{  sum— subset! j J ;  subset! j 1-0;  j+-l;  } 


-7 


- 

/*  updates  points  and  their  values  */ 
void  update (void) 

<  int  i.  j; 

if  (newval  <  value[l]  ) 

{  for  (i-0;  1<5;  1++) 

for  {j-0;  j<ns;  j++)  resultl [i] C j J-cresult [ i ] [ j J ; 
if  {gap— 1) 

{  value [2 l«value[ 1 1 ;  value [ 1 l*newval ; 
pntt2i-pnt[ll 5  pnt[ 1 1-newpt j 

} 

if  (gap— 2) 

{  value [0 ]«value [ 1 ] i  value [ 1 ]-newval ; 
pnt [0i-pnt[ 1 ] ;  pnt ( 1 ]-newpt; 

} 

} 

else 

{  if  (gap— 1) 

{  value[0]-newval;  pnt(0]“newpt;  } 
if  (gap— 2) 

{  value  1 2 ]*newval ;  pnt ( 2 ] -newpt;  } 

} 


- 

/*  does  a  binary  enumeration 
void  binen(lnt  p) 

{  int  XX,  j; 

xx-p;. 

for  (j"0;  j<ns;  j++)  subset[ j  1-0*, 
j-0; 

while  (  XX  >  1  &&  j  <  ns  ) 

{  subset[j]  -  xx%2: 

XX  -  (xx-(xx%2) )/2;  j+-l; 

} 

subset [ j ]-l : 


- 

/*  restarts  the  line  search  Iterations  */ 
void  restart ( void) 

{ 

int  1,  j,  k,  p,  pi,  temp,  flag,  sum; 

if  ( value( 1 ]<tempv)  /*  save  current  best  value  */ 

{  tempv-value [ 1 ] : 
for  (i-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  aresult(l] ( j ]-resultl[i] [ j ] ; 

} 

P-0; 

if  (jj--ns-l)  /•  if  full  set  is  to  picked  up  next  •/ 

{  for  (j-0;  j<ns;  j++) 

{  subsetIjJ-1;  p-(2*p)  +  1  ;  } 
ptvalue ( ) ; 
if  (newval  <  tempv) 

{  tempv  -  newval: 
for  (i-0;  i<5;  i++ ) 

for  (j-0;  j<ns;  >++)  aresult[ 1 J [ j ]-cresult [1 J [ j ] ; 

} 

if  (jj<ns-l) 


{  P-0; 

for  {  j-0;  j<jj+l;  j++)  t*  left-point  and  corresponding  subset*/ 

.  subset[j]-l;  p-2*p+l;  } 

lor  (  j-jj+i;  j<ns;  j++)  subset[jl-0; 
ptvalue();  valuel01-newval;  pntl0]-p;  tempv-newval; 
for  {i*0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  aresult [ i J [ j J-cresult [ i 1 [ j ] ; 
for  (  j-0;  j< (ns- ( j j+1 ) ) ;  j++)  subset[j]-0;  /•  right-point  and  */ 
for  (  j-(ns- ( j j+l) ) ;  j<ns;  j++)  subset! j]-l;  /*corresponding  subset*/ 
for  (  j-jj+1;  j<ns;  j++)  p  *-  2; 
pnt [ 2 ] -p ;  ptvalue ( ) ;  value [ 2 ] -newval j 
if  (newval<teiapv) 

{  tempv-newval; 

for  (i-0;  1<5;  i++) 

^  for  (j-0;  j<ns;  j++)  aresultli] t j l-cresult[l] [ j ] ; 

pnt[l]»(pnt[0]+pnt[2] )/2;  binen(pntt 1 ] ) ;  sum-0; 
for  (j-0;  j<ns;  j++)  sum+-subset [ j 1 ;  j-0; 

while ( sum< j  j  +  1 ) 

{  sum+-l-subset[ j ] ;  subsetrjl-i;  j+-l;  } 
while(sum>jj+l) 

{  sum--subsetl j ] ;  subset! j]-0;  j+-l;  } 
ptvalue();  value ! 1 ] -newval ; 
for  (i-0;  i<5;  i++) 

for  (j-0;  j<ns;  j++)  resultl ! i ] ! j J-cresult ( i ] ! j ] ; 


fn!2J-'S';  fn!3J-'T';  > 
fnI2]-'M';  fn!31-'B';  } 


. . / 

/*  prints  the  output  on  the  screen  */ 
void  output(void) 

{  int  i.  j; 
char  fn!4]; 
if  (objflag— 1) 

{  fn!0]-'C';  fn!l]-'o' 
if  (objflag--2) 

{  fn!0]-'T';  fn!l]-'l' 

If  (  tempv  <  value! 1]  ) 

{  value ! 1 ]-tempv; 
for  (1-0;  i<5;  i++) 

^  for  (j-0;  j<ns;  j++)  resultl!!]! jj-aresult!lll j J; 

printf("  OUTPUT  STATEMEHTXn" ) ; 

- - - NnXnXn"); 

printf(  Input  FILE'S  name  is  '\s’  \n\n\n",  fname); 

PrlntfC  Corresponding  to  \c\c\c%c  value  of  %5.2f,  the  f ollowing\n" . 

IJ*  fn!2I,  fn!3],  value!!]); 

table  gives  details  of  recommended  surveillance  sensorNn"); 
prlntf(  mechanism  deployments\n\n" ) ; 

Ptintf(- _ 

printf ( "Sensor\tRadius\t\tLocation\t\tCost\t\tTime\n\n" ) ; 

for  (j-0;  j<ns;  j++) 

{  if  (  resultl!41 ! j 1  >  -i  ) 

prlntf(-  %d\t%5.2f\t  (»5.2f.  %5.2f.  %5.2f)\t  %4.2f\t  %5.2f\n' 

J+1.  sensor!0J!j],  resultl!0] ! j ] ,  resultl ! 1 ]! j ] .  resultl (2 ]! j ] . 

^  resultl!3]! j],  resulol !4] ! j ] ) ; 


> 

/* 


/*  checks  whether  the  give'  space  may  be  covered  by  a  given 


*/ 

*/ 


/*  set  of  sensor  mechanisms »  and.  If  so.  computes  the  'value'  *t 

/*  of  the  cover  relative  to  given  objective  function  */ 

void  select (void) i 

void  sort(vold): 

void  xyz(vold); 

void  check(vold)i 

Int  p.  quit; 

Int  pointer [LIM] ,  slct[LIM].  part[LIM]; 
float  X.  y.  z.  costlmeiLIM] ; 
void  ptvalue( void) 

{  Int  j.  count;  /*  selects  deployment  */ 

count-0; 

while  (count  <  CHANCES) 

{  for  (j-0;  j<ns;  j++) 

{  cresult[0] [ j ]-0;  cresult[l] [ j ]-0;  cresult[2] [ j ]-0; 
cresultis] [ j i-0;  cresulti4i [ j i»-2;  part[j]-0; 

} 

select();  sort();  p-1;  qult-0;  newval-0; 
while  (  p  <  pleces+1  &&  quit  --  0  ) 

{  xyz();  check(); 

If  (quit  --  0)  P+-1; 

} 

If  (quit  --  0)  count-CHANCES+1;  /*  a  cover  found  */ 
else  count+-l; 

} 

If  (qult--l)  /•  not  a  cover  */ 

{  newval-BIG; 

for  (j-0;  j<ns;  j++) 

{  cresultI0] I j ]-0;  cresult 1 1 1 [ j 1-0;  cresultI2 ] i j 1-0; 
cresult[3 j [ji“0;  cresultI4jljj--2; 

> 

} 


- - 

/*  selects  a  deployment  of  a  given  set  of  sensor  mechanisms  */ 
void  select(vold) 

{  Int  j,  r.  nu«; 

for  (j-0;  j<ns;  j++) 

{  num-seed*7;  r-sensor [ 1 ] [ j  ] ; 

slct[ j  ^-( subset[ j 1 ) * (num%r) ;  seed+-l ; 

} 

} 


- 

/*  sorts  deployed  sensor  Tiechanlsms  by  costs  or  by  time  */ 
void  sort(vold) 

{  Int  1.  j.  k,  col.  tempi,  count; 
float  temp; 

If  (objflag--!)  k-3;  else  k-4; 

/•  copying  Into  costlme  */  1-0; 
for  (j-0;  j<ns;  j++) 

{  If  (subset ( j I--1 ) 

{  col-lndex[  j  I-t-slct  [  j  1 ; 

costlme [lj-deploy[ki [col] ; 
polnterilj-j ;  1+-1; 

) 

/*  sorting  */  count-jj; 
while  (count>0) 


{  for  (i"Cs  i<count;  i+-+) 

{  if  (costlme[l]>costlme[l-«-ll ) 

{  teinp«costlme [  1  ] :  costlme [ i ]>coBtlme [  i-i-1  ] ; 
costlme[ 1+11-temp;  templ-polnter[i] ; 
pointer[i ]-pointer[ i+1 1 ;  pointer! 1+1 ]-  tempi; 

> 

) 

count--l ; 

} 

} 

- 

f*  computes  the  origin  of  the  piece  corresponding  to  p  */ 
void  xyz(vold) 

{ 

Int  num.  a,  b.  c,  temp; 
num-p-1;  temp  -  l*w; 
c  -  num/temp;  num  %-  temp; 
b  -  num/l;  a  -  num%l: 

temp-space [ 3 ] :  x-a*temp;  y-b*temp;  z-c*temp; 

) 


/ . . */ 

/*  verifies  that  the  point  (x,y,z)  is  covered  */ 
void  check(vold) 

{  float  a.  b.  c,  r.  sa.  sb.  sc.  sr; 
int  i,  j,  col,  flag; 
fl.ag-0;  i-0; 
while  (flag»-0  &&  i<jj) 

{  j-pointer[il ;  col-index[j]  +  slct[jl; 
r-sensor[0] [ j ] ;  a-deploy[0] [col ] ; 
b-deploy( 1 ] [col ] ;  c-deploy[2 ] [col] ; 

sr-r*r;  sa-(a-x)* (a-x) ;  sb-(b-y)*(b-y) ;  sc- ( c-z ) * ( c-z ) 
if  {(sa+sb+sc)  <-  sr) 

{  flag-l;  cresult[0] [ j ]-a;  cresult[ 1 ] [ j ]-b; 

cresult [ 2 ] [ j ] -c;  cresult [3] [ 1 J -deploy [ 3] [col ] ; 
cresult[4I [ j i-deploy[4] [col] ; 
if  (part[jl— 0) 

{  if  (objflag--!)  newval+-cresult[3 J [ ) ] ; 
if  (objflag--2  &&  newval<cresult [4 ] [ j ] ) 
newval-cresult[4] [ j ] ; 
part[ j  J-1; 

} 

> 

i+-l ; 

} 

if  (flag--0)  quit-1; 


} 
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